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I ' Abstract 

^: 

5^ . We consider different extensions of the standard model which can 

give rise to the small active neutrino masses through seesaw mecha- 
nisms, and their mixing. These tiny neutrino masses are generated 
^ . at some high energy scale by the heavy seesaw fields which then get 

OO I sequentially decoupled to give an effective dimension-5 operator. The 

renormalization group evolution of the masses and the mixing param- 
eters of the three active neutrinos in the high energy as well as the 
ly-^ ■ low energy effective theory is reviewed in this article. 
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1 Introduction 



1.1 Neutrino oscillations: The current status 



The field of neutrino pliysics lias made immense progress in the last decade, 
which was initiated when the Super-Kamiokande (SK) experiment in Japan 
[1] reported the evidence for oscillations in the atmospheric neutrinos. Now 
there is compelling evidence that solar, atmospheric, accelerator and reactor 
neutrinos oscillate, which implies that the neutrinos are massive and the 
leptons mix among themselves. 

The atmospheric neutrinos are produced in the Earth's atmosphere by 
cosmic rays. The flux of cosmic rays that lead to neutrinos with energies 
above a few GeV is isotropic. Hence one expects the downward and the 
upward-going fluxes of multi-GeV neutrinos of a given flavor to be equal. 
The underground SK detector found that for multi-GeV atmospheric muon 
neutrinos the zenith-angle dependence deviates from this expectation and 
the deviation can be explained when one invokes z/^ — ?• v^ oscillations. The 
oscillations of muon neutrinos into other flavors have also been confirmed by 
the energy spectrum obtained from the controlled source experiments K2K 
[2] and MINOS [3J. The allowed region for the oscillation parameters, Am^^j^^ 
and sin^26'atm, is shown in Fig 1(a) 
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Figure 1: (a) The region of the atmospheric oscillation parameters Am^^j^j 
and sin^ 2^atm obtained from the SK, K2K and MINOS experiments [4j; 
(b) The allowed region in the neutrino oscillation parameter space from solar 
neutrino data and KaniLAND experiment ^. 



As can be seen from the figure, the MINOS data is capable of measuring 
^^atm with high precision, while SK put stronger bound on sin^ 26'atm- The 
results from the short-baseline (SBL) experiments (like CDHS [6], NOMAD 
[7] etc.) show that the u^^ — >■ Ug oscillations can be present only as small 
sub-dominant effects and also put strong bounds on the active-sterile mixing 
angles in z/^ — )■ Ug oscillations, an oscillation channel whose sub-dominant 
effect is not yet ruled out completely. 

The pioneering solar neutrino experiment by Davis and collaborators us- 
ing ^^Cl reported a solar electron neutrino flux significantly smaller than 
that predicted by the standard solar model, and this deficit in the number of 
electron neutrinos is known as the "solar neutrino problem" . The puzzle per- 
sisted in the literature for about 30 years, and then the charged current (CC) 
and the neutral current (NC) data from the SNO experiment [S], combined 
with the SK solar neutrino data [9], provided direct evidence for neutrino 
oscillations in solar neutrinos. However, four different solutions were there 
to explain the solar neutrino oscillations [10]: (i) the LMA or large mixing 
angle solution (Am^ = 5.0 x 10"^ eV^, tan^ 6*0 = 0.42), (ii) the low mass 
solution (Ami, = 7.9 x 10~^ eV^, tan^ ^0 = 0.61), (iii) the vacuum solution 
(Am| = 4.6 X 10"^° eV^ tan^ 6^0 = 1.8) and (iv) the SMA or small mixing 
angle solution (Am| = 5.0 x 10"® eV^, tan^ 6q = 1.5.10"^). The results from 
the controlled source experiment KamLAND [11] confirmed the LMA solu- 



tion and ruled out the other three possibilities. Fig 1(b) shows the allowed 
region of the solar neutrino oscillation parameters Am0 and tan^ ^0. 

Combining the results obtained from the solar, atmospheric and the reac- 
tor neutrino oscillation experiments described above, the current knowledge 
about the neutrinos is that there are three neutrino flavors (z/^, a G {e, /i, r}) 
which mix to form three neutrino mass eigenstates (z/,,? G {1,2,3}). These 
mass eigenstates are separated by Am^^- = mj — m| where, rriij denote mass 
eigenvalues with i,j G {1,2,3}. The two sets of eigenstates are connected 
through z/q = {UpMNs)ai^i, where t/pMNS is the Pontecorvo-Maki-Nakagawa- 
Sakata neutrino mixing matrix [121 1131 El US] in the basis where the charged 
lepton mass matrix is diagonal. This mixing matrix is parametrized as 

f/pMNS = i^-W-Q, (1) 

where 

W = ?723(^23,0)f/i3(^i3,^)f/i2(^i2,0), Q = Diag{e"^'^S e-^<^^ 1} . (2) 

Here Uij{6, 6) is the complex rotation matrix in the i-j plane, 6 is the Dirac 
CP violating phase, (pi are the Majorana phases, and P is the flavor phase 
matrix (Sometimes the flavor phases are called as the unphysical phases since 





Best fit 


3cr range 


Am^i 10-^eV^ 


7.65 


7.05 - 8.34 


Aml^ [lO-^eV^ 


2.40 


2.07- 2.75 


sin^ 6^12 


0.304 


0.25-0.37 


sin^ ^23 


0.50 


0.36 - 0.67 


sin^ ^13 


0.01 


< 0.056 



Table 1: The present best-fit values and 3cr ranges of oscillation parameters 
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they do not play any role in the phenomenology of neutrino mixing or beta- 
decay.) Finally, with all the above definitions, U takes the form 

/ C12C13 S12C13 sise"*'' 

U = -C23S12 - S23Si3Ci2e*'' C23C12 - S23Sl3Si2e*'' S23C13 1 , (3) 



C12C13 


S12C13 


sise 


-C23S12 - S23Sl3Ci2e*'' 


C23C12 - S23Sl3Si2e*'' 


S23C13 


S23S12 - C23Sl3Ci2e*'^ 


-S23C12 - C23Si3Si2e*'^ 


C23C13 



where Qj and Sij are the cosines and sines respectively of the mixing angle 6ij . 
The current best-fit values and 3a ranges of these parameters are summarized 
in Table [H It is still not known whether the neutrino mass ordering is 
normal (mi < m2 < m^) or inverted {m^ < mi < 1112). Many other high 
precision oscillation experiments are going on and also being planned in order 
to measure the neutrino oscillation parameters with higher accuracy and to 
determine the neutrino mass ordering. 

As can be seen from the PMNS parametrization of the neutrino mixing 
matrix in Eq. Q, the angle 9i^ plays a crucial role in the determination of the 
Dirac CP phase 6. As shown in the Table [H ^13 can also be consistent with 
zero at 3cr. However, this data also implies that assuming the error to scale 
linearly upto 3cr within the physical range of sin^ ^13, there is a hint of ^13 > 
at ~0.9(T. It has been shown that the solar and KaniLAND data implies a 
non-zero ^13 at ~1.5(j [13, [H]. But when combined with atmospheric, long- 
baseline reactor and CHOOZ data, the significance is lowered since the hint 
for a non-zero ^13 from the atmospheric data is not so robust and depends 
on the details of event rate calculations and the treatment of theoretical 
uncertainties IT8l. 



1.2 Absolute masses of the active neutrinos 

While the neutrino oscillation experiments are not sensitive to the abso- 
lute neutrino masses, the beta decay and the neutrinoless double beta decay 
(Oi//3/3) processes are. At the same time, it is possible to estimate ^^ rrii 



from cosmology also. In case of beta decay, the non-zero neutrino mass 
would modify the Kurie plot, regardless of whether the neutrinos are Dirac 
or Majorana particles. The effect will depend on m^ = (^^ \Uei\^mf) , and 
if the neutrino masses are small, it will be visible only near the end point 
of the Kurie plot. The Mainz [20j experiment has placed the upper limit of 
mp < 2.3 eV. The upcoming beta-decay experiments like KATRIN [21] will 
be sensitive to uip > 0.2 eV and will thus improve the bound by an order of 
magnitude. The Oi//3/3 decay, on the other hand, is sensitive to the effective 
Majorana mass of the electron neutrinos, defined as mee = \ YliUei''^i\j ^^^ 
will be observed only if the neutrinos are Majorana particles. A non-zero 
signal for the Oz//3/3 decay will put bound on the specific combination of the 
neutrino masses and the Majorana phases given by rriee- The current limit 
put by the Heidelberg-Moscow experiment [22j is rrige ^ 0.9 eV. The cosmic 
microwave background radiation (CMBR) carries the imprint of the neutrino 
masses since in the standard Big Bang model, for the standard model (SM) 
interactions of the neutrinos, the neutrinos are abundant like the photons till 
the epoch of nucleosynthesis when they decouple from the thermal bath of 
the photons. It is also possible to get information about the neutrino masses 
from the study of the large scale structure as an active neutrino species of 
mass nil, will tend to wash out all structures upto a scale ~ l/m,^ by free- 
streaming. Recent results from the Wilkinson Microwave Anisotropy Probe 
(WMAP) and the surveys on the large scale structure put the limit ^2- irti < 
0.67- 2.0 eV [221 IM|. 

The very fact that the active neutrinos are massive demands an exten- 
sion of the SM. In the framework of the SM, since there is no right-handed 
neutrino, the neutrinos are massless at the tree-level, and they cannot have 
a Dirac mass even at loop level. So the only other possibility is the lepton 
number violating Majorana mass term. But lepton number is a symmetry of 
the SM, though accidental, and if that symmetry is to be obeyed, Majorana 
masses also cannot be generated at loop level. It can also be seen that the 
Planck scale (Mpi) effect cannot introduce the required neutrino mass in the 
SM as it can only generate a neutrino mass ~ C(w|;^/Mpi) ^ (9(10~^eV), 
and hence cannot explain the atmospheric mass squared difference. Hence 
generally the neutrino masses are incorporated at the tree-level by adding 
new fields to the SM at high energy scales. The most favored mechanisms 
to generate such small neutrino masses are the so called seesaw mechanisms 
which need the introduction of one or more heavy fields, while maintains the 
SU(3)cX SU(2)iX U(l)y gauge group structure of the SM. There are also 
other models like Inverse seesaw [2S] , the model with a singly charged singlet 
proposed originally by Zee [26l[27|, the model with a doubly charged singlet 



^27\ f2E\ . etc. Recently another new model has been proposed in [2S], where 
a pair of vector like leptons and also a Higgs quadruplet are added to the 
SM to generate neutrino mass. However, some models, like the Zee's model, 
cannot predict neutrino mixing parameters consistent with the current data. 
We will discuss some of these models of neutrino masses in detail in Section [2j 

1.3 RG evolution of neutrino parameters 

Since the neutrino mass is generated at the high scale while the neutrino 
masses and mixing parameters are measured experimentally at a low scale, 
the renormalization group (RG) evolution effects need to be included. The 
current experimental data in Table [T] shows that in the neutrino sector two 
of the three mixing angles are large, while the third one is small, which is 
rather different from the quark sector where all three mixing angles are small. 
Because of the large values of the two mixing angles, RG evolution of the 
neutrino masses and the mixing parameters plays an important role in the 
neutrino sector, which is not the case with the quark sector. RG evolution 
will be even larger if the neutrinos happen to be quasi-degenerate. 

The radiative corrections in different theories of neutrino masses are ex- 
pected to be different since the heavy particles couple differently to the SM 
fields present. However below the mass scale of the lightest of the heavy par- 
ticles the effect of all heavy degrees of freedom are integrated out to get an 
effective theory of neutrino masses. The RG evolution of neutrino masses and 
mixing parameters in the low energy effective theory as well as in different 
high energy theories will be discussed in Section |3] and Section H] respectively. 

2 Generation of light neutrino masses 

2.1 Low energy effective theory of neutrino masses 

The low energy effective Lagrangian needed to explain the non-zero active 
neutrino masses can in general be expressed as a series of non-renormalizable 
operators, the dominant one being the dimension-5 operator given as [30] 



-Ck ~ K^bhlLftxp ■ (4) 

where II and (j) are respectively the lepton and Higgs doublets belonging to 
the SM. Here ^5 is the effective coupling which can be expressed in terms of 
a dimensionless coupling 05 as ^5 = as /A with A some high energy scale. In 
this picture the SM serves as an effective theory valid upto the mass scale 
A, which can be taken to be the mass of the lightest of the heavy fields. 
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Figure 2: Generation of the Majorana neutrino mass from the low energy 
effective Lagrangian given in Eq. (jl]) after spontaneous symmetry breaking. 



However, the specific form of K5 will depend on the high energy field content 
and the interactions present at the high scale. The operator shown in Eq. (jl]) 
violates lepton number by two units and gives rise to Majorana masses for 
neutrinos after spontaneous symmetry breaking, laii, ~ ^k^v"^, as shown in 
the Fig. 121 Here v is the vacuum expectation value (vev) of the Higgs field 
such that 

Spontaneous / U 



Taking v 
05 ~ 1. 



symmetry breaking \ /^ 

246 GeV, a neutrino mass of ~ 0.05 eV implies A 



(5) 
10^5 GeV if 



2.2 High energy theories: Seesaw mechanisms 

There are four possible ways to form a dimension-5 gauge singlet term as 
given in Eq. (jl]) at low energy through the tree-level exchange of a heavy 
particle at the high energy: (i) each li-(j) pair forms a fermion singlet, (ii) 
each of the Il-^l and 0-0 pair forms a scalar triplet, (iii) each [^-(j) pair forms a 
fermion triplet, and (iv) each of the Il-Il and 0-0 pair forms a scalar singlet. 
Case (i) can arise from the tree-level exchange of a right handed fermion 
singlet and this corresponds to the Type-I seesaw mechanism [311 ESI [331 
[311 [35] ■ Case (ii) arises when the heavy particle is a Higgs triplet giving 
rise to the Tj^e-H seesaw mechanism [271 [3B] ■ For case (iii) the exchanged 
particle should be a right-handed fermion triplet, which corresponds to the 
Type-Ill seesaw mechanism [39l[l0j. The last scenario gives terms only of the 
form 1^2 ei, which cannot generate a neutrino mass. We describe the three 
different seesaw mechanisms in Section I2.2.1f[^?^751 in detail. A summary of 
the form of k and hence the effective light neutrino mass at the low scale in 
different types of seesaw is given at the and of this section in Table [2 There 
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is another model, similar to the seesaw models, that can predict the light 
active neutrino masses and known as the inverse seesaw model. This model 
will be discussed in Sec I2.2.4[ for the sake of completeness. 

2.2.1 Type-I seesaw 

The simplest extension of the SM to incorporate small active neutrino mass 
is to introduce right-handed singlet fermions Nji in the theory, which are 
singlets under the SM gauge group. Hence these Nji fields are essentially 
right-handed neutrinos. Presence of these new fields allows new terms in the 
Lagrangian 

Cm = ^NmN - biU^N - (nY^^Ul + h.c.) , (6) 

where II andj/) are respectively the lepton and Higgs doublets belonging to 
the SM and = ia'^(j)*, a"^ being the second Pauli matrix. Here we do not 
write the generation or the SU(2)£, indices explicitly. The field N is defined 
as A^ = Nfi + N^, where A'"^ is the CP conjugate of the right-handed field 
Nfi. Yn is the Yukawa coupling for the singlet fermion and Mn is the mass 
matrix. Thus the complete Lagrangian of the theory becomes 

C = CsM + Cn , (7) 

and after spontaneous symmetry breaking it is possible to write the neutrino 
mass terms as 

-C...4(prW)(j:j; m:)(;^)+-, (8) 

where 

mz, = {v/y/2)Y^ (9) 

is the Dirac mass matrix for the neutrinos generated after the electroweak 
symmetry breaking when the Higgs gets the vev f , as given in Eq. ([5]). Thus 
the complete mass matrix for the neutrinos becomes 

which when block-diagonalized gives the eigenvalues (see Appendix for deriva- 
tion) 

mi ^ -mz)MN"^mJ , (11) 

ma ^ Mn , (12) 

9 





The effective vertex 
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Table 2: Summary of tlie low energy effective couplings and the effective 
neutrino mass in,^ = — ^k in the three seesaw scenarios. Here, 1n(^s) 
are the Yukawa couplings for the heavy singlet (triplet) fermion present in 
Type-I(Type-III) seesaw and Mn(Ms) is the mass matrix (A^ = Nji + N^; 
S = S/j + S^). In Type-II seesaw, Ma is the mass of the heavy triplet 
Higgs, Ya is its Yukawa coupling with the SM lepton doublet l^, and Ag is 
its coupling with the SM Higgs (p. 



where we have assumed that Mn ^ m^), i. e. the eigenvalues of Mn are much 
larger than the eigenvalues of m^ and kept terms upto 0(jnD/M-^). Thus 
Eqs. (TT2|) and (TTT]) show respectively that eigenvalues of the matrix 102 are 
large, while those of mi are small and hence the eigenstates corresponding 
to these small eigenvalues should serve the purpose of the mass eigenstates 
of the light active neutrinos. Thus the presence of the heavy right-handed 
neutrinos will produce the light active neutrino masses and this mechanism 
of making one particle light at the expense of making another one heavy is 
called the seesaw mechanism. The seesaw obtained by adding these heavy 
right-handed singlet fermions to the SM is called the Type-I seesaw. 
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In the low energy limit we have an effective theory described by 

l^n = «:/,(/£V£0)(0Vdi)+h.c., (13) 

= -l^fg yLc^Jib4>dj - {SacSbd + SabScd) + h.C. , (14) 

where k is a symmetric complex matrix with mass dimension (—1) and e = 
ia"^ is the completely anti-symmetric tensor in the SU(2)2, space. Generation 
indices f,g ^ {1,2,3} are shown explicitly and a,b,c,d G {1,2} are the 
SU(2)i indices. In writing Eq. ( IT^ we have used 

{<y^)ab{<^')cd = '2'SadSbc - SabScd 
=^ {cr''£)baicr'e)dc = ^Eda^bc - £ba£dc (15) 

and utilizing the (pd ^ 4>a symmetry, we can write 

'^^da^bc — ^ba^dc = " {^ab^dc + ^db^ac) ■ (16) 

The relevant diagrams in the complete theory giving rise to the effective 
operators in the low energy limit are shown in the topmost row in the Table [21 
The "shaded box" on the left hand side of the equivalence in the middle 
column represents the effective low energy vertex k, while A(^a) and A(^b) are 
the amplitudes of the diagrams labeled as (a) and (6) on the right hand side. 
The amplitudes are given by 

Aia) = lfi'{Y^M^^Y^)^^eca£bdPL, (17) 

A(^b) = ifJ.'{Y^M^^Y^)^^e,debaPL, (18) 

with e = 4 — D where D is the dimensionality that we introduce in order to 
use dimensional regularization. Note that A(^b) is obtained from A(^a) just by 
d ■v^ a interchange. Using Eq. flTSjl one finally gets 

A(a) + A^b) = -«/U' (Fi^Mn^Fn)^^ {SabScd + SacSbd) Pl ■ (19) 

This is equal to the left hand side of the figure mentioned with the identifi- 
cation 

K = 2Y^M^'Y^ , (20) 

as shown in the Table [2J From Eqs. (1141) and (1201) . one gets the neutrino 
mass after spontaneous symmetry breaking to be 

m, = -^Y^M^'Y^ (21) 

which is the Type-I seesaw relation. As the energy changes, the heavy sin- 
glets get decoupled one by one at their respective mass scales and start 
contributing to the light neutrino mass through the effective operator. 

11 



2.2.2 Type-II seesaw 

In the Type-II seesaw, we consider the SM extended by a charged Higgs 
triplet transforming in the adjoint representation of SU(2)l 



A 






A+/v^ A++ 
AO -A+/V2 



(22) 



where A++ = (A^ - i/\^)/V^, /SP = (A^ + iA'^)/V2 and a' = {(j\(r^,a^ 
are the Pauli matrices. Following the notation of [HI |12], the Lagrangian is 
given by 

^ = -^SM + ^A , (23) 

where 



Here 



'-'Ajkin 



^A — ^A,kin + ^A,(l> + ^A,Yukawa ■ 



Tr 



{D^Ayo^A 



(24) 



(25) 



Ai 
2 



-M^Tr (At A) - -^ [Tr (A+A)] 

A2 

' 2 



[Tr ( A"^ A) ] - Tr ( A"^ A A"^ A) - A40V Tr (A"^ A) 



-A5,^'[A',A] 



7!^ 



«cr2 



+ h.c. 



(26) 



C 



A,Yukawa 



-^{Y^)f,fJC{^a,)Aii + h.c., 



(27) 



where C is the charge conjugation matrix with respect to the Lorentz group. 
The covariant derivative of the Higgs triplet is given by u 



(2^ 



Di^A = d^A + t\l-giB^A + tg2 [W^, A] , 



where gi and g2 are the U(l)y and SU(2)i gauge couphngs respectively. 
With the interactions shown in Eqs. (I24j) - 027j) . after electroweak symmetry 
breaking the triplet Higgs A will get a vev given by (Aq) ~ Agf ^/2v2M^. 
This triplet Higgs vev will also contribute to the gauge boson masses and 



^ We use GUT charge normalization: | (g^- 



3 /'„GUT^ 

5 



5i 



SM^ 
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will alter the p-parameters from the SM prediction p ~ 1, at tree level and 
hence will get a strong constraint from the current precision data [44j- 

Once the triplet Higgs A gets its vev after spontaneous symmetry break- 
ing, the Lagrangian in Eq. (127|) produces the neutrino mass term as 

C-A,Yukawa = ^Ia^^l (Ao)z/L + h.C. , (29) 

and thus using the expression for (Aq), the neutrino mass is given as 






™- = —^7^ ■ (30) 



^A 



In case of Type-II seesaw, only one diagram in the complete high energy 
theory contributes to the effective low energy neutrino mass operator, as 
shown in Table [2], and we have 

•^ = ^^nT(^^)fa(^<^<:^bd + £ab£cd) , (31) 

and comparison with Eq. f ll4p gives 

2A6>A 



Mi 



(32) 



Hence finally one gets the neutrino mass to be 



li 



THy = 1^ = ^— 1 (33) 

4 2 M2 ' ^ ' 



which is the same as the Type-II seesaw relation, as given in Eq. (15U]) . 

Just like the right-handed neutrinos in case of Type-I seesaw, the Higgs 
triplets in Type-II seesaw will decouple step by step at their respective mass 
scales and the effective theories have to be matched against each other. The 
decoupling of the right-handed neutrinos only contributes to the effective 5- 
dimensional neutrino mass operator, while the decoupling of the Higgs triplet 
also gives a contribution to the SM model Higgs self-coupling because there 
is a coupling between the SM Higgs doublet and the Higgs triplet given in 
Eq. fl26|) . The matching condition for the Higgs self-coupling at the threshold 
is given as 



lAe'^ 

M2 



A'^'^^ = A + 2^. (34) 
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2.2.3 Type-Ill seesaw 

Type-Ill seesaw mechanism is mediated by heavy fermion triplets transform- 
ing in the adjoint representation of SI]{2)l and has been considered earher 
in [391 SO] • Very recently there has been a renewed interest in these type of 
models. The smallness of neutrino masses usually implies the mass of the 
heavy particle to be high ~ 10^^~^^ GeV, as shown in Chapter 12.11 How- 
ever, it is also possible that one or more of the triplets have masses near 
the TeV scale, making it possible to search for their signatures at the LHC 
|15| HHl SZl HH So]. In such models, the Yukawa couplings need to be small to 
suppress the neutrino mass, if no fine tuning of the parameters is assumed. 
Lepton flavor violating decays in the context of Type-Ill seesaw models have 
also been considered in [50] . Recently it has also been suggested that the 
neutral member of the triplet can serve as the dark matter and can be in- 
strumental in generating small neutrino mass radiatively [5T] . 

In the Type-Ill seesaw, there are right handed fermionic triplet S/j added 
to the SM at the high scale which is singlet under U(l)y, while transform 
as a triplet in the adjoint representation of 811(2)^. This triplet can be 
represented as 



'R 



^VV2 s+ \ _ J:],a 



s^, -j:%/V2 J v^ 



(35) 



where S^ = (S)j =i= iS|.)-\/2- For the sake of simplicity of further calculations, 
we combine S/j with its CP conjugate S^ to construct 

S = S^ + Sg . (36) 

Clearly, E also transforms in the adjoint representation of SU(2)j^. Note that 
though formally S = E*-^, the individual elements of E are not all Majorana 
particles. While the diagonal elements of E are indeed Majorana spinors 
which represent the neutral component of E, the off-diagonal elements are 
charged Dirac spinors. 

Introduction of this triplet field will introduce new terms in the La- 
grangian. The net Lagrangian is 

C = CsM + C^ , (37) 

where 

£s = ^T.,kin + ^J:,mass + ^T.,Yukawa ■ (38) 
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Here, 

£s,fcm = Tr[a|E] , (39) 

Cj^^mass = -^TrpMsS], (40) 

C-J,,Yukau,a = -IZv^F^^S^ - ^^^^Ev^Fs^L • (41) 

Here we have not written the generation indices exphcitly. Mg is the Majo- 
rana mass matrix of the heavy fermion triplets and Yy, is the Yukawa cou- 
phng. Since the fermion triplet S is in the adjoint representation of SU(2)i, 
the covariant derivative of S is defined as 

D^T. = d^T. + ig2[W^,Y], (42) 

where g2 is the SU(2)l gauge coupling. Unlike A, S being a singlet under 
U(l)y does not couple to i?^. 

The new term £s in the Lagrangian, as shown in Eq. (1551) . can be ex- 
panded as 



+92 (w^T^YPr"^ + W^T^^Pl^ + h.c.) - g2W^W^ 
_^Ms^ - QsjMsSOp + h.c. j 

- (0°s5Fs^^l + V2<J)'^Yj,Il + 0+s5Fs/l - y20+^yj^ + h.c. 

(43) 

Here we have defined the four component Dirac spinor 

vl/ ^ E+^ + S- , (44) 

for our convenience, while the neutral component of Sij is still in the two 
component notation. In Eq. f H3|) . the first two lines come from C-£,kin, the 
third line corresponds to the Majorana mass term in C-^^mass and the terms 
in the last line corresponds to the Yukawa coupling terms in CY,,Yukawa, as 
given in Eqs. (139|) - (HT|) . After the electroweak symmetry breaking, the mass 
matrix for the neutral fields become 

where mo = (v/y/2)Y^ is the Dirac mass matrix of the neutral fields. Thus 
the mass matrix in Eq. fHSjl looks the same as that obtained in Eq. flTOj) 
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and hence for large Ms, diagonalization of the mass matrix will produce 
light active neutrino states via seesaw mechanism, as obtained in Sec I2.2.1I 
The seesaw achieved here with the help of the neutral component of the 
fermionic triplet is known as the Type-Ill seesaw mechanism. Eq. (H5l) also 
implies that there will be a mixing between the light and the heavy neutral 
states, however the mixing angle will be (9(]m£)/Ms) and hence very small 
for large Ms- 

Since the heavy fermion triplets added to the SM at the high scale have 
charged components also, they will modify the masses of the charged leptons 
belonging to the SM, in addition to the generation of the small active neutrino 
masses. With the addition of the triplet fields, the mass term of the charged 
lepton sector after electroweak symmetry breaking becomes 

= - (i^ *;)>'=( ^'J-iu *:)>';( i; 



R 
R 



(46) 



where iul is the Dirac mass matrix of the SM charged leptons and 



denotes the complete mass matrix for the charged leptons. Eq. fH7|) shows 
that the inclusion of the charged fermions as components of the heavy triplets 
does not change the masses of the charged leptons of the SM upto the order 
0((]m£), ]mL)/Ms)- However, there will be mixing between the states Il-^l 
and Ir-^r, but the mixing angle is small in the ioid, ihl <^ Ms limit. The 
correction to the charged lepton masses due to the charged components of 
the triplet fermion in C([(iin£), tul) /M.^]'^) can be calculated easily from [17] . 
As can be seen from the Table [21 the diagrams in the complete Type-Ill 
seesaw theory giving rise to the effective operators in the low energy limit 
are very similar to the case of Type-I seesaw. Here the amplitudes Ai^a) and 
A{h) are given by 

A(a) = ^/i^(lfMs'>s)^J(£V)„fe(£V)e,]PL, (48) 

^(fe) = ^/x^(lfM^^Fs)^J(£V),5(5V),JPi. (49) 

Using Eq. (TT5|) one finally gets 

A^a) + A(p) = -ifi' (Y^M^^Ys) ^^ {eab^cd + £ac£bd) Pl , (50) 
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which gives 

K = 2yjMs^Fs • (51) 

From Eqs. JHM and (I5T1) one gets the neutrino mass after spontaneous sym- 
metry breaking to be 

Mi,= -^YiM^'Yj, (52) 

which is the Type-Ill seesaw relation. Here, v denotes the vacuum expecta- 
tion value of the Higgs field. 

2.2.4 Inverse seesaw 

Apart from the three types of seesaws described in Sec. 12.2. 111^X51 there is 
another well-known scenario known as the Inverse Seesaw. In the Inverse 
Seesaw scenario [25j, additional SM gauge singlets are introduced together 
with a small Majorana mass insertion through the additional right handed 
heavy singlets which explicitly breaks the lepton number. The minimal ver- 
sion of this Inverse Seesaw scenario requires the addition of two right-handed 
neutrinos iV^ and two left-handed SM gauge singlets S^, where / G {1, 2} is 
the generation index [551 ES]- We construct the field A^-^ and S-^ as 

Nf = nI^ + N^^ , (53) 

S^ = Si + S^^, (54) 

and then the Lagrangian for the theory will be given by 

^ = -^SM + ^is ) (55) 

where Cjs is the contribution from the fields added to the SM to produce 
the Inverse Seesaw and is given as 



^IS — ^IS,kin + ^IS,Yukawa + ^IS,mass , (56) 



where 



^isMn = I'N^mf.N'^ + l'Sfmf.S^, (57) 

ClSXukan^a = -N^ (Y^) jg^^li + h.C. , (58) 

jO.IS,mass = -S\Mr) ^^N^ -^'S^flfgS''^ + h.c. , (59) 

where the generation indices g, f are written down explicitly. Here /x is a 
complex symmetric 2x2 matrix and In and M.^ are arbitrary 3x2 and 2x2 
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matrices, respectively. Without loss of generality, one can always redefine the 
extra singlet fields and work in a basis where n is real and diagonal. After 
the electroweak symmetry breaking, the Lagrangian giving rise to the 7x7 
mass matrix for the neutral fields u^, N-^ and S-^ can be expressed as 

/OmoOWz/fX 

C3(l7ZN^S2)\^l M^ iV^ +h.c. , (60) 

^ ^ \ Mn fi J \S^, J 

where the Dirac mass matrix is defined as mid = {v/V2)Y^. It should be 
noted that M/j and fi being the mass terms of the SM singlet fields do not 
depend on the scale of the SU(2)i symmetry breaking. 

At the leading order in m/^M^^, the active light neutrino mass matrix is 
given by 



m,. 



mDM^V(M5) V5 = F/iF^, (61) 



where F = hibM^^. In this case, for fi ~ 10^ eV, the light neutrino mass 
can be 101,^ ~ 0.01 eV if F ~ 0.3 x 10^^. Thus for the Yukawa couplings 
In ~ 0.1 - 1, the heavy singlets can be in the mass range 10^-10^ GeV and 
the seesaw scales can be lowered by orders of magnitudes compared to the 
Type-I seesaw. 

3 RG evolution of neutrino masses and mix- 
ing in effective theories 

Now we consider the radiative corrections to the masses, mixing parameters 
and couplings in the effective low energy theory of neutrino masses. As can 
be understood from the discussions in Sec. [2|, there is a unique dimension- 
5 operator given in Eq. (j4j), that gives rise to the light neutrino masses 
after spontaneous symmetry breaking and hence is the same for all three 
types of seesaws. Thus the RG evolution equations will depend solely on the 
underlying theory. Throughout this paper, we consider the SM as the low 
energy effective theory. However, to discuss the RG evolution in the effective 
theory, we will also consider the Minimal Supersymmetric Standard Model 
(MSSM) and discuss the effect of tan/3. 

In order to evaluate the /3-functions, one need to compute the renormal- 
ization constants for wavefunctions, couplings, etc. For this purpose any 
regularization scheme can be chosen from dimensional regularization, Pauli- 
Villars method. Ultra- Violet cutoff etc. and then the scheme for renormaliza- 
tion has to be fixed [60]. For the purpose of calculation of the renormalization 
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Yu[ZY^Yu + YlYd + Ti 



(2Ti- 
sfl^i 

13 „2 
15^1 



5 6*1 



Qgl)K 
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l,9l - 3gl - fgi 
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-3) 



Tr 



^^ u ^ u 



Tr 



y'ye + syjyd 



Table 3: Evolution equations for the Yukawa couplings Fe, Yu^ Yd, gauge 
couplings (71, 5^2, 5'3, Higgs self-coupling A (SM only) and the effective neutrino 
mass operator n in the SM and the MSSM [361 EH EHl 1104] . Here we have used 
the GUT charge renormalization and hence gi = (jfUmfied _ ^5^3 ^^sm^ This 
convention has been followed through out this review. Here Px = n{dX/dfi). 



constants, the gauge also has to be fixed. However, the final /3-functions must 
be independent of the particular regularization as well as the renormalization 
scheme used for the calculations, and also of the gauge choice. /3-functions 
can be determined from the renormalization constants using the functional 
differentiation method, as described in [361 EO]- The running equations for 
the Yukawa couplings, gauge couplings, Higgs self-coupling (in case of the 
SM only) and the effective neutrino mass operator in the SM and the MSSM 
are given in Table [31 Finally, after electroweak symmetry breaking, the light 
neutrino mass matrix is given by 



m. 



v^ 
T'^' 



(62) 
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and thus will have the same evolution as n. To simplify later discussions, the 
evolution equation for ni;^ can be expressed as 

167r^/3iii^ = P^Tm„ + ui„P + a„m.^ (63) 

where Px = dX/dln{fi/GeV) and 

P = CeYX- (64) 

Here the values of Ce and a^ depend on the underlying theory, and can be 
read off from Table [3] when the theory is the SM or the MSSM. Without any 
loss of generality we can always choose the charged lepton Yukawa matrix 
Yf. as well as the quark Yukawa matrices Yu and Y^ to be diagonal at the 
high scale. Then from the RG equations in Table [3] we get that they will 
remain diagonal at all energy scales, and so will P. Thus the evolution of 
the components of k,, and hence of in,^, will be proportional to themselves. 

3.1 Evolution equations for neutrino parameters from 
matrix equations 

At any energy scale /i, the neutrino mass matrix nij^ and the charged lepton 
Yukawa Y^Y^ can be diagonalized by unitary transformations via |lU4j 

t/j,(Ai)^mj.(/i)f/^(/i) = Diag(mi(/i),m2(/i),m3(/i)) , (65) 

[/e(/i)ty;y,(/i)[/,(/i) = Bmgiy'M,yM,y'M) , m 

where Ui, and Ue are unitary matrices and the neutrino mixing matrix will 
then be given by 

f/pMNs(/i) = Ulifi)UM . (67) 

Since in the effective theory Y^Y^. remains diagonal at all energies, as dis- 
cussed in the previous section, ?7e(yu) = 1 in Eq. f pB]) and from Eq. fIBTl) one 
has L^pmns(a*) = Uy{ij). Thus the RG evolution of f/pMNs(A*) will be governed 
by the running of nij, only. The evolution of the mixing matrix Uv will be 
given by [61] 

dU,, 

-^^U.T, (68) 

where t = ln(/i/GeV)/167r^ and T is an anti-Hermitian matrix defined as [M] 

167r2lmr,, = -!I^i^^ Imi^' (70) 

mj + rrij ■' 
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where P' = UlPU^, with P defined by Eqs. flB5]) - flM|) . In order to ob- 
tain the RG evolution equations for the mixing angles and the phases, one 
has to solve the system of nine coupled equations in the parameters C,k = 
{Oi2, O23, Oi3, 6, (pi, 4)2, Se, Sfj,, Sr} , obtained from Eq. (168|) using the definition 
of T as given in Eqs. (l69l) - (170l) and then using the parametrization of f/pMNS 
as given in Eqs. ([I])-©. 

3.1.1 RG evolution of the mixing angles and phases 

The RG evolution equations for the mixing angles can be written in general 
as [HmfTMIHT] 

^"^ +Ax + 0(^13) , (71) 



X 



9 



13 



where X G {Ou, 023,613, 6, 4>i,(p2}- The differentiation is performed w.r.t. 
t = ln(yu/GeV)/167r^. It can be seen that the quantities Dx = for all X 
except D13. Evolution of the mixing angles are given by 



A12 
A23 
Ar3 



Cyl .^^^ 2 I "^le^**^^ + m2e 



2i<t>2 



2 

Cyl 



sin26'i2S 



23 



Am^ 



sin 29' 



23 



77126 



"12 






+ m3 



Am 



atm 





(72) 


mie^*'^! + 7773 ^ 


,(73) 


A"^atm(l + C) 



sin 26*12 sin26'' 



7713 



X 



Cyl 

—— bm zt7i2 bm z,U23 -r- 2 / 1 AX 

2 AmL„,(l + C) 

[7711 cos (201 — (5) — (1 + C)m2 cos (202 — 5) — C"^3] 



(74) 



where C = '^"^©/^^tm ^^"^ % ~ cos 9ij, Sij = sin 9ij. As can be seen from 
Tabled C = Ce = -3/2 for the SM and C = Ce = 1 for the MSSM. Here y^ 
and 7/^ are neglected compared to y"^. The quantities governing the evolution 
of the Dirac CP phase 6 are given as 



As 



2C7/2|^l!!j4sin(20i-202) 



7773 



Amf,^(l + C) 



[c23 (miCi2 sin {25 - 20i) + 7712(1 + C)-Si2 sin {25 - 202)) 



D, 



+ cos 26*23 (777iSi2sin20i + 7772(1 + C)ci2sin202)] 



■ sin 2^12 sin 29 



^3 
''A7nL^(l + C) 



[ttt-i sin (201 — 5) — {1 + C)^2 sin (202 — 5) + (m^ sin 5] 



(75) 



(76) 
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while those for the Majorana phases (pi, 02 are 

2/ ofl "^iS?2sm20i + (l + C)m2C?2sin202 
A^^ = 2Cy^'lm3 cos 2623 . 2 ^ , /-^ 

mim2C^2^23 sin (20i - 202) ) 

Ami /' ^''^ 



^. mis|2sm20i + (l + C)m2C?2sin: 

As^ = 2Cy^< 1713 008 2023 



Amf,^(l + C) 

mim2sl2sl3 sin (20i - 202) ) ,^^. 

Am| J ■ 

3.1.2 RG evolution of the light neutrino masses 

Using the definition of U^ in Eq. f l65|) and T in Eq. fl68|) . the RG evolution 
of the light neutrino masses is obtained from Eq. f lB^ to be 

rhi = (Rea^ + 2Re/^J m^ , (79) 

where no summation over the repeated index 'i' is to be taken. The evolution 
of the individual masses in terms of the mixing parameters becomes 

mi = [a^ + Cyl{2sl^sl3 + Gi)]mi, (80) 

rn2 = [a, + Cyl{2cl2sl3 + G2)]m2, (81) 

7713 = [ay + 2Cylcl3cl3] ms , (82) 

where 

Gi = —Si3 sin 2^12 sin 26*23 cos 5 + 2S13C12C23 , (83) 

G2 = Sia sin 26*12 sin 26*23 cos 6 + 2S13S12C23 . (84) 

From the Eqs. fl5U]) - fl5^ it can be seen that the evolution of a particular 
mass eigenvalue is proportional to itself upto 0(6*13) and thus if some rrii is 
zero to start with along with 6*13 = 0, it will remain so. However, it is the 
characteristic of the 1-loop RG evolution only, and breaks down when the 
2- loop contributions are taken into account [S^. 

To study the RG evolution of the neutrino masses and the mixing param- 
eters in the effective theory, we consider /^o to be the high energy scale below 
which the effective theory gives the correct description of the light neutrino 
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masses, which we also take to be the mass of the hghtest heavy particle re- 
sponsible for the seesaw mechanism. Then at any energy scale /i, the value 
of the mixing angles can be expressed as 

% = el+ [ AJ{t')dt' + 0{9^3) (85) 

J to 

^ ^0 +%A. + 0(A.^i3,A2), (86) 

where to = ln(/io/GeV)/167r^ and 6*°- is the value of the angle at the high 
energy /xq- In Eq. (186|) . A,- is defined as 

in the SM, where g2 is the S\J{2)l gauge coupling, whereas rrir and Mw are 
the r lepton and W boson masses respectively. In the MSSM, 

1 / \ 2 / 

MSSM _ J- / fl'2"^r \ f^ , ^„_2 a\ 1„ / /^O 



Numerically, one has A™ ^ — 1.4 x 10"^ when fiQ = 10^^ GeV and /i = 10^ 
GeV. For MSSM, A^^sm ~ _i.3 x 10-^(1 + tan2 /3), where jj = 10^ GeV and 
tan /3 can take values upto ~ 50. Hence in both the cases one can treat these 
quantities as small parameters. The quantities kij can then be written from 
Eqs. (^^-(^M as 



C I mie^*"^! + mae^*'^^ 

-Sm 2^12523 —^ 

2 Am, 



ku = --sin 2^124^ ^^ ^, (89) 

"^0 



C 

^23 = -—sin 26*23 



mae^*''^^ + m3 I „ I rriie^^'^^ + ms 



2- 



-^12 A 2 ~^ "^12 



Amftm Amf,„Xl + C) 



(90) 



=13 = —sin 26*12 sin 26123-: — 2"^^ TV ^ 

2 Am2^^(l + C) 



[mi cos (201 - 5) - (1 + C)m2 cos (202 - 5) - C^^s] • (91) 

The same results are obtained in [631 EH [65] following a slightly different 
approach given in [BSl EZl ES]- Similar integrated evolution equations can be 
written for the Majorana phases as 

0i = 0° + A;^,A. + O(A,6^i3,A^), (92) 

where 0° is the value of 0j at yUo and k^^ can be read off directly from Eqs. fl77|) - 
rHj) . However, the running of the Dirac CP phase 5 has to be considered 
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carefully, since D^ is non-zero and 6i^ is allowed to take small values including 
zero. This issue will be discussed in detail in Sec I3.1.3I 

From ki2 in Eq. (189|) it can be seen that the solar mixing angle 9i2 gener- 
ically has the strongest RG effects among the mixing angles. The reason for 



this is the smallness of the A/tIq associated with it, in particular compared to 



^'^atm) which leads to an enhanced running for quasi-degenerate neutrinos 
and for the case of an inverted mass hierarchy. The running is maximum for 
\(t>i — </'2| = 0, and minimum for |0i — 02| = 7r/2. As it is clear from Eq. fl89l) . 
the direction of the running depends solely on CA,-. Hence for evolution from 
a high to a low energy scale, 9i2 always increases in the MSSM, and decreases 
in the SM. From Eq. fl90|) it is evident that the direction of the ^23 evolution 
depends on CA,- as well as on the hierarchy. However, the running of ^13 
depends on specific combinations of the CP phases, as shown in Eq. (1911) . If 
the symmetry ^13 = is implemented at the high scale /zq [6^, which is the 
case for many neutrino mass models, the maximum ^13 value can be achieved 
with the choice 

20i-<5o = O, |202-'5o|=vr, (93) 

5q being the Dirac CP phase at /xq and finally 

^max < l^sin2^i2sin2^23^r^[mi + (l + C)m2 + |C|m3] . (94) 

To consider the running of the Majorana phases, one gets combining 
Eqs. dUD-dZH]) and Eq. (El 

kd,^-ks2 = 2C cos 2^12823— — 2~ sin (2(/)i - 2(/)2) , (95) 

AttIq 

which shows that if (0i — ^2) = at some scale, it will remain so at all energy 
scales, upto 0{6i^). Moreover, if (0i — ^2) is small at some scale so that we 
can write sin (201 — 202) ~ 2(0i — 02), the running of (0i — ^2) is proportional 
to itself. However, the 0(^13) term may become important for large tan/3 
values in case of the MSSM and then it will be possible to generate (0i — ^2) 
radiatively 

The running of the mass eigenvalues is significant even in the SM or for 
strongly hierarchical neutrino masses due to the factor a^ in the RG evolution 
equations given in Eqs. (I80l) - (l82l) . As can be seen explicitly, the evolutions 
are not directly dependent on the Majorana phases, and the dependence on 
the Dirac CP phase 6 is proportional to sin0i3. Moreover, apart from the 
MSSM with very large tan/3 or at very high energy values, the running of 
the mass eigenvalues is solely controlled by the term proportional to a„, and 
there will be very small dependence on the mixing parameters. Thus in such 
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cases, the running is given by a common scaling of the mass eigenvalues 
and can be given by 



rrii ^ m^ exp 



au{t')dt' 

to 



(96) 



Some generic features of the RG evolution of the light neutrino masses and 
the mixing parameters in the effective theory have been studied extensively 
in literature [^11 ESI EHl EHl ESI E] • These effects can have interesting conse- 
quences such as the generation of large mixing angles [721 [731 EH EHl UM EI] > 
small mass splittings for degenerate neutrinos [78l[79ll8nl[8n[82ll83l[8il[85l 
[861 187] , or radiative generation of ^13 starting from a zero value at the high 
scale [521 [SSI [SSI [SO] • Some specific features of the RG evolution, like the sta- 
bihty of mixing angles and masses [SSI [SD [S21 [SSI [SI] , possible occurrence of 
fixed points [951 [SSI [SI] have also been studied. RG induced deviations from 
various high scale symmetries like tri-bimaximal mixing scenario [Ml l65l 198] 
or quark- lepton complementarity [SSI [SHI 11001 HOll I102J and correlations with 
low scale observables have also been explored in detail. 

3.1.3 A subtlety at ^13 = 

As mentioned already, Eq. f[751) clearly suggests that As and hence 6 diverges 
for ^13 — )• 0. This problem is overcome by requiring that D^ = at 6'i3 = 0, 
which gives the following condition on 6 at 6*13 = [6T] : 



micos20i - (1 + C)m2cos202 -C"^3 /^i^x 

cot = ■ -. — ■ . (97) 

nil sm 201 — (1 + C)^^2 sm 202 

The above prescription works for the calculation of evolution when one starts 
with vanishing 613. However on the face of it, it seems to imply that the CP 
phase 6, which does not have any physical meaning at the point 613 = 0, 
should attain a particular value depending on the masses and Majorana 
phases, as given in Eq. f[^ . Moreover, getting the required value of 6 pre- 
cisely when ^13 = would seem to need fine tuning when one starts from 
some non-zero ^13, unless this value of 5 is a natural limit of the RG evolution 
when 6*13 — )• 0. The problem also propagates to the evolution of 613, since 
y4i3 in Eq. (17^ depends in turn on 6. The evolution of all the other param- 
eters, viz. rrii, 6^12, 6*23 and 0j is independent of 6 upto 0{9i^), and hence 
will have continuous, non-singular evolution even at ^13 = 0. This apparent 
singularity in 6 has been explored in [62j by analyzing the evolution of the 
complex quantity Wes, which stays continuous throughout the RG evolution 
and shows that a fine tuning is indeed required, but that is to ensure that 

25 



6'i3 exactly vanishes. However, if the parameters happen to be tuned such 
that ^13 vanishes exactly, then the limiting value of 6 as ^13 — t- is always 
the one given by the prescription mentioned in Eq. (1971) . 

Even after understanding the origin of the apparent singularity in the 
evolution of S and hence of ^13, a necessity still remains to have a clear 
evolution of parameters that reflect the continuous nature of the evolution 
of elements of the neutrino mixing matrix t/pMNS- This can be achieved by 
choosing the basis as Pj = {rrii, 612, 6*23, Of^, (pi, Jqp, Jcp} where the quantities 
JcPyJcp are defined as 



-^CP = 2'^12Cl2S23C23-5l3C?3Sin(5 , (98) 



•^CP = o'5l2Cl2S23C23Sl3q3COs5 , (99) 



instead of the conventional basis Vs = {rrii, 612, 623, 613, 0j, 5}. From Eqs 
and dnnD it is seen that Jcp, J'cp — )■ as 6^13 — )■ and thus are well-defined. 
The RG evolution equations for Jcp and J^p are given as 

JcP = Aj + 0(613), (100) 

j^p = A'j + o{ei3), (101) 



with 



m3 



mi sin 201 — (1 + C)^^2 sin 202 



Aj — C'Z/r ■512^12-523 '^23 7 2 /-, , an ' (102) 



A"^ftm(l+C) 



1713 



nil cos 201 — (1 + C)^^2 cos 202 — C'^3 

^J = Cy^Si2Ci2S23C23 T 2 /-, , ^n .(103) 

^"^atm(l + C) 

In the new basis Vj, the RG evolution for 6*13 is considered instead of 
^13, as is traditionally done. This quantity turns out to have a nonsingular 
behavior at ^13 = 0. Moreover, since ^13 > by convention, the complete 
information about 6*13 lies within ^13. Also, the possible "sign problem'^ of 
6^13 is avoided. In terms of the new parameters Jcp and J'q^, the RG evolution 



^Usually the convention used in defining the elements of L/pMNS is to take the angles dij 
to lie in the first quadrant. Ues can then take both positive or negative values depending 
on the choice of the CP phase S. In the formulation of Eq. (l74|) the sign of A13 can be such 
that 9i3 can assume negative values during the course of evolution and in such situations 
one will have to talk about the evolution of 16*13!. Our formulation in terms of 6*^3, as 
shown in Eq. (|104p . naturally avoids this problem. 
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equations for Of^ becomes 

Oh = All + 0{el,), (104) 

All = SC'l/rT — 2 /, , an I-^cp [mi sin2(/)i - (1 + C)'^2 sin2(/)2] 

+Jcp [rrii cos 20i — (1 + C)"^2 cos 202 — C^^s] r • (105) 

Thus the evolution equations in basis Vj are all non-singular and continuous 
at every point. In particular, even when 6 shows a discontinuity, Jqp as well 
as Jqp change in a continuous manner. This very fact can be used to write 
down the approximated integrated evolution equations for Jqp, Jqp as 

Jcp = j0p + fcj^pA, + O(A,^i3,A2), (106) 

^cp = Jc°p + fcjtpA. + 0(A,^i3,A2), (107) 

where J^p, Jqp are the initial values at /io and kj^^, kji can be obtained 
from Eqs. (ll02p - (11031) . From the Jqp, Jqp values, the Dirac CP phase 5 can 
be determined unambiguously at any energy scale. 

4 RG evolution of neutrino masses and mix- 
ing in high energy seesaw models 

As discussed in Sec [31 the RG evolution of the neutrino masses and mixing 
parameters in the low energy effective theory is the same for all three types 
of seesaw scenarios and depends only on whether the low energy effective 
theory is the SM or the MSSM. However, this is not true in the high energy 
theory, when the heavy particles responsible for seesaw remain coupled to 
the theory. Hence in this case the RG evolution of the different neutrino 
parameters should also depend on the interaction of these heavy particles 
with other fields. The importance of including the effects from energy ranges 
above and between these mass thresholds when analyzing RG effects in GUT 
models has been pointed out in [721 EDI ED EDI [IHSl EZl 1101 [Ml [IDS] • These 
effects are typically at least as important as the RG evolution effects from 
below the thresholds since the relevant couplings may also be of order one. 

The diagrams contributing to the renormalization constants of the dif- 
ferent quantities at high energy are shown in [36], [12] and [107] in case of 
Type-I, Type-II and Type-Ill seesaw respectiveljQ. Finally the RG evolution 



^ [m 132] actually considered Type-I + Type-II case where they have one heavy right- 
handed fermion and a triplet Higgs added to the SM (and the MSSM also). However, we 
will consider the three types of seesaw scenarios separately. 
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equations for the charged lepton Yukawa matrix Y^. and the heavy particle 
Yukawa matrix Yx and (here X = N for Type-I, X = A for Type-II and 
X = E for Type-Ill) can be given as 

167r2/3y^ = Y,F + a,Y,, (108) 

2^ J YxG + axYx Type-I & III 



where F and G are defined as 

F = D,YX + DxYJ,Yx, (110) 

a = B,YX + BxY^Yx. (Ill) 

The quantities De,Dx, Bf.,Bx and Og, «x with the SM as the low energy 
effective theory are given in Table HI The /9-function for the heavy particle 
mass Mx is given as 



167r2/3M^ = I Mx {YxY^x) + {YxY^x) ^x + «MxMx Type-I & III , 

(112) 



«Mx ^x + «M^ Mx Type-II , 



with a-Mx &^d a'jy^ defined by 

r Type-I , 

«M, = < 4Ai+A2 + Tr[r2rA]-|^?-6(7| Type-II, (113) 

I -12^2 Type-Ill, 

a'^^ = 2A4m2 + -|A6|2 (Type-II), (114) 

where m is the bare mass of the SM Higgs and AjS are the couplings 
associated with the triplet Higgs A, as given in Eq. f l26|) . 

After electroweak symmetry breaking, these coupled heavy fields will con- 
tribute to the generation of the light neutrino mass via seesaw, as given in 
Sec l2.2.Hl2l2l3] for the three different seesaw scenarios 

r -^Fn'^MnVn Type-I , 
^^=\ +T^ Type-II, (115) 

[ --^Y^M-^^Yj, Type-Ill . 

In this section we would like to study the radiative corrections to the quantity 

Y^M^^Y^ Type-I , 
Q^{ -^ Type-II, (116) 

F^Ms^Fs Type-Ill , 
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which is a well-defined quantity at different energy scales to denote the con- 
tribution to the light neutrino mass matrix from the coupled heavy fields and 

2 

finally gives the light neutrino mass matrix as jm.,^ = —\Q after spontaneous 
symmetry breaking. 

In the complete theory when the heavy fields are coupled, the RG evolu- 
tion of Q can be obtained from the running of Yx and Mx (running of Ag 
is also required in case of Type-II seesaw, which can be obtained in [HI US]) 
and finally the evolution of Q can be written as 

167r% = QPQ + P^Q + aQPQ, (117) 

with Pq defined as 

Pq = C,Y}Y, + C'xY},Yx. (118) 

The quantities Cg, C'x and aq for three seesaw scenarios with the SM as the 
low energy effective theory are also given in Table HI Finally m^ can be 

2 

obtained from Q as im,y = —\Q, after spontaneous symmetry breaking. 

In case of Type-I seesaw, the heavy fermion singlets do not have any 
gauge interactions and hence the presence of these singlets does not affect the 
running of the gauge couplings. However, for the other two seesaw scenarios, 
the RG evolution of the gauge couplings gi, g2, gs will depend on the number 
of heavy fields present. Let n be the number of heavy fields present at some 
energy scale. Then the evolution of the gauge couplings can be written as 

167^2/3,^=M^ (119) 

where the values for 6jS in the three seesaw scenarios are tabulated in Table [51 
As given in the Table none of the heavy fields has any strong interactions 
and hence 63 is always the same as its SM value. The Higgs triplets present 
in Type-II seesaw have Y = 1 and thus couple to both U(l)y and SU(2)i 
gauge fields, while the triplet fermions in Type-Ill seesaw have Y = and 
thus couple only to the SU(2)l gauge fields. 

We do not give the evolution equations for the different Higgs couplings 
or for the up- and down-type quark Yukawa couplings here. The evolution 
equations for these couplings can be obtained in [611 SD 1107] for the different 
seesaw scenarios. 

4.1 Sequential decoupling of heavy fields 

The most general case of the high energy theories would have been the one 
when there are any arbitrary number of right-handed singlets in Type-I see- 
saw, or any number of triplet scalars in case of Type-II seesaw or arbitrary 
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De 


Dx 


Oe 


Type-I 

Type-II 

Type-Ill 


3/2 
3/2 
3/2 


-3/2 

3/2 

15/2 


Ti - U - U 

Tn - U - U 
Tin - U - U 




Be 


Bx 


ax 


Type-I 

Type-II 

Type-Ill 


-3/2 
1/2 
5/2 


3/2 
3/2 
5/2 


^I ~ 2q5i - 4^2 
T' -^^2 9J2 

Tin - M - f 5i 




C'e 


C'x 


ag 


Type-I 

Type-II 

Type-Ill 


-3/2 
1/2 
5/2 


1/2 
3/2 
3/2 


^rj. y 2 9 2 

T,, - 2T[, - Zgl + A + /(Ai) 
2^/// - ^A - hi 




Ce 


Cx 


Oi>^ 


Type-I 

Type-II 

Type-Ill 


-3/2 
-3/2 
-3/2 


1/2 
3/2 
3/2 


2Ti + \- Sg'i 
2Tn + A - 3gl 
2Tiii + A - 3gl 



Table 4: The quantities defining the running of Ye, Yx, Q and n for three 
types of seesaw scenarios [611 HH 1107] . Here T/ = Tr[y'j"y'e + ^n^n + 
3Y^Y^ + SYiYd], Tn = Tr^tF, + 3i;ty„ + 3^]^^, T;, = Triy^yA] and 
Tin = Tr[i;tFe + ^Y^Y^ + SY^Y^ + SFJfJ- The function /(A^) is defined as 
/(A-) = -8Ai - 2A2 - 4A4 + 8A5 - (4A4m2 + | Aep) MX^ where m is the bare 
mass of the SM Higgs doublet and Aj are the couplings associated with the 
triplet Higgs A, as defined in Eq. (12^ . 





h 


b2 


63 


SM 

Type-I 
Type-II 
Type-Ill 


41/10 

41/10 

41/10 + 3n/5 

41/10 


-19/6 

-19/6 
-19/6 + 2n/3 
-19/6 + 4n/3 


-7 
-7 
-7 

-7 



Table 5: 6, in three types of seesaw scenarios with the SM as the low energy 
theory [SH 1121 1107] . Here n is the number of heavy particles coupled to the 
theory at any particular energy scale. 

number of fermion triplets in Type-III seesaw and these heavy fields decouple 
one by one at different thresholds. 

Let us first consider the most general case of Type-I and Type-III seesaw 
when there are r heavy fields (singlets in case of Type-I and triplets in case 
of Type-III) having masses Mi < M2 < ■ ■ ■ < Mr—i < Mr. We consider 
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a quantity R which contains the contribution from the coupled as well as 
decoupled heavy fields at any energy scale and is also well-defined at all n, 
and then finally gives the light neutrino mass matrix as 



m. 



-R 



(120) 



after electroweak symmetry breaking. 

Above the heaviest mass Mr, all the r- fields are coupled to the theory 
and contribute to R as 



(r+l) (r + 1) 

R = Q, 



(121) 



(r + l) 



where Q denotes the contribution from the r coupled heavy fields and is given 

by 



(r+l) (r + l) (r + l) (r+l) 

g = F| M^i Yx (/i > Mr) , 



(122) 



(r + l) 



where X = N for Type-I and X = E for Type-Ill seesaw. Here Yx is a 
[r X Up] dimensional matrix {np is the number of flavors, which is 3 in our 
case) given as 



(r+l) 

Yx 



\ iyx)r,i 



iyx)i,np \ 

iyx)r,nF J 



(123) 



(r + l) (r+l) (r + l) 

M.X is a [r X r] matrix and Q as well as R is a [np x np] dimensional matrix. 
We use the super-indices just to keep track of the number of coupled flelds. 
Below the scale M^ , the heaviest of the heavy fields decouples from the theory. 

(r) 

Integrating out this degree of freedom gives rise to an effective operator k. 
The matching condition at /i = Mr is 



(r) 



Mr 



(r+l) 

2(>^J) 



{Mr 



(r+l) 

' (yx) 



rj 



(124) 



Mr 



where no summation over 'r' is implied and i,j G {1,2, •• ■ jUp}- This con- 
dition ensures the continuity of R at /i = Mr. In order to get the value of the 
threshold Mr, we need to write the above matching condition in the basis 
where Mx = Diag(Mi, M2, ■ ■ ■ , Mr). Here it is worth mentioning that the 
matching scale has to be found carefully since Mx itself runs with the energy 
scale, i.e. Mi = Mi{n). The threshold scale Mi is therefore to be understood 
as Mi(/i = Mi). 
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In the energy range M^-i < /x < M^,, R will be given as 



(r) 

IR 



L{r) 



(r) 



^^ + Q. 



(125) 



The first term in Eq. (11251) is the contribution of the integrated out heavy 



(r) 



fermion of mass Mr through the effective operator k. The second term rep- 
resents the contribution of the remaining (r — 1) heavy fermions, which are 

(r) (r) 

still coupled to the theory. Mx is now a [(r — 1) x (r — 1)] matrix while Yx 
is a [(r — 1) X np\ dimensional matrix given as 



Yy^ 



iyx)r-i,i 



iyx)i,np 



iyx)r-l,n. 



\ 



\ 











/ 



(r) 

Yx 



Mj. integrated out . 



(126) 



(r) (r) (r) 

Finally Mx and Yx constitute Q, which is [np x np] dimensional. The 
matching condition at yU = Mr_i is 



(r-l) 

Kij 



(r) 



Mr 



Mr 



+ 2{Yx\r-i) (Mr-i)-' [Y: 



X)l^r-l)j 



,(127) 



Mr^ 



where no summation over '(r — 1)' is to be taken. 

Generalizing the above sequence, we can say that if we consider the in- 
termediate energy region between the {n — 1)**^ and the v}^ threshold, i.e. 
Mn > n> Mn-i, then all the heavy fields from masses Mr down to Mn have 

(n) 

been decoupled. In this region the Yukawa matrix Yx will be [{n — 1) x np] 
dimensional that couples the (n — 1) coupled fields with np flavors and will 
be given as 



Yx^ 



( {yx\i ■ 


• (2/^)l,n^ 


{yx)n-l^ ■ 


■ (l/x)„_i,„p 








^ 






(n) 

= Yx, 

heavy fermions with 
masses M„ — Mr 
integrated out . 



:i28) 



(n) 



Mx will be [{n — 1) x (n — 1)] dimensional matrix involving the mass terms 
of all the coupled heavy fields. In this energy range IR will be 









(129) 
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EFTl - 



R=i« 



EFT2 



,,, (2) (2) 

K, Yx. Mx 



R = i^ + Q 

(2) (21 (2) (21 



EFT3- 



K, Yy, Mx 



R^ 



1(3) 



* — "FULL" theoiy 
Yx, Mx 

R = Q = Y^Mx^Yx 



(3) (31 (31 (31 

Q = Y^Mx^Yx 




Figure 3: Sequential decoupling of the heavy fermions and construction of 
R at different energy scales, for r = 3. Here X = N in case of Type-I and 
X = E for Type-Ill seesaw. Finally the light neutrino mass matrix is given 

2 

as nijy = —^11 after spontaneous symmetry breaking, where v is the vacuum 
expectation value of the SM Higgs. 



with 



(n) (n) 



(n) 

M 



X 



(n) 

Yx 



(130) 



Note that R, k and Q are [np x np] matrices. The matching condition at 

/x = Mn is given by Eq. fll27p with r replaced by (n + 1). 

At low energies fi < Mi, when all the heavy fields are decoupled, Q{fi) = 
(1) (1) (1) 

Q{n) = and ]R(/i) = (l/2)K(/i). Fig. |3] shows the expressions for R at 

different energy scales for the case of three heavy fermions i.e. for r = 3. 

2 

Finally the light neutrino mass matrix 101,^ is obtained as ni;^ = —'y^, after 
the electroweak symmetry breaking. 

In case of Type-I and Type-Ill seesaw, the concept of sequential decou- 
phng is important since atleast two heavy fields are needed to generate the 
observed pattern of the light neutrino mass matrix and the heavy fermions 
can have non-degenerate masses in general. However, the case of Type-II 
seesaw is much simpler since only one heavy scalar triplet is sufficient to 
give rise to the small masses of the three active neutrinos and their mixings. 
Hence in this case one has only one threshold at fj, = Ma and R at the two 



33 



different energy regimes wiii be given by 

R=|? i"*^;-)- (131) 

\ i/t (/i < Ma) , ^ ' 

wtiere Q and k are defined in Eq. flll6p and Table [2] respectively. Tlie matcli- 
ing conditions at yU = Ma will be given by 



«:.. L. = -^^A6(1a)„- , (132) 

Ma 



2 

IMa ~ ^' IMa "^ j^ l''-^! I Ma 



^ Ima - ^ Ima + ¥72-|^6P L. • (133) 



Finally the light neutrino mass matrix nit, is obtained as tsiIu = —-^11, after 
the electroweak symmetry breaking. 

For the sake of convenience and to be consistent with the existing litera- 

2 

ture, we will refer to in,^ = —y^ ^s the effective light neutrino mass matrix 
at any energy scale, in the rest of the thesis. However, it must be understood 
that the quantity 'f ' is present only after electroweak symmetry breaking 
and hence, strictly speaking, this relation is valid only in that energy regime, 
while R is a well-defined quantity at all energy scales. 

Finally, the RG evolution of the light neutrino mass matrix will be given 
by 

167r^/3m, = -P'^m^ + m^P + a^-m^, , (134) 

which is of the same form as Eq. (1631) . Here P = Pq, a^ = aq for /i > M3 
and P = Pk, Ciu = O-k. for yU < Mi. In case of Type-I and Type-HI seesaw, the 
running of the light neutrino masses in between the thresholds will be given 
by the running of both Q and k, as shown in Eq. (11291) . For the energy scale 

Mn > /i > M„_i, the relevant quantities will be Q and k. The RG evolution 

(n) (n) 

of Q can be obtained from Table H] with the substitution Yx -^ Yx and 

(n) 

Mx — )■ Mx- In case of Type-II seesaw there is only one threshold and the 
running of all relevant quantities contributing to the running of the neutrino 
mass matrix can be read off from Table HI 

It can be seen from the interaction of the heavy fields that they do not 

(n) 

contribute to the 1-loop correction of the effective operator k, even when 

(n) 

they are coupled to the theory. So the evolution of the effective vertex k can 
be written as 

IOtt^A") = PjK+^K'P. + aJK^ (135) 
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where 

P. = CeY^Y. + CxY^Yx, (136) 

and the quantities Cg, Cx and ««; can be obtained from Table H] and are 
independent of ra, the number of heavy fields coupled at any energy scale. 

4.2 RG evolution of neutrino mixing parameters 

In order to evaluate the RG evolution of the light neutrino masses and the 
mixing parameters from the matrix evolutions, we proceed in the same way 
as done in Sec 13.11 Without the lose of generality, we choose to work in the 

basis in which M^ and YJY^. are diagonal at the high energy. 

As stated in the last section, for yU > M3 and yU < Mi, the evolutions of 
mij, will be given by the evolution oiQ [P = Pq) and k, {P = P^) respectively 
and hence the evolution of the angles, phases and light neutrino masses can 
be given in simple analytic forms. P and F, in Eqs. fll34p and fllOSp . are 3x3 
matrices with the rows and columns representing generations. We denote the 
elements of P and F by Pfg and Ffg. If we write the evolution equations in 
the basis Vs, the apparent singularity at ^13 — )■ will be present, as can be 
seen from [l2l[57j. As already discussed in Sec 13.1.31 this singularity can be 
removed using the basis Vj [107| . Hence we discuss the RG evolution of the 
mixing angles, phases and the light neutrino masses in the Vj basis in the 
following sections. 

4.2.1 Evolution of mixing angles 

Running of the two large mixing angles 612 and 623 in the basis Vj, as given 
in Table is also the same as that in the Vs basis since the quantities Sij 
and Qfj, defined as 

^13- Am^,^(l+C) ' ^23- A^^;^ '^12" A^^^ ' U^O 

Q mim:j sin2(/)i <-• r?i2r»3 sin2(/)2 c m,im2 sin {2<f>\ — 2ij>2) /-. qoN 

•^IS - Aml,^{l+0 ' '^23 - Aml,^ ' ^^12- A^^^ ' (138) 

depend on the mass eigenvalues and Majorana phases only, and not on the 
Dirac CP phase 5. However the running of 9\.^, as seen from the Table [6], 
depends on the quantities Af,, Bf, defined as 



-_i_ 4 {ml + m\) J^p ± 8mim3( J^p cos 20i + Jcp sin 20i) 

a^^ltm (1 + C) 



>^i3 ~ „A^2 n , A\ ' U'jyj 



j± _ 4 {ml + ml) Jgp ± 8m2ms{J^p cos 202 + Jcp sin 2^2) , , 

35 





327r2^i2 


327r2 ^23 


Pu 


Q+ sin 2^12 





-'22 


-Q|2sin2^i2C23 


(2^3C?2+2|3^?2) sin 2^23 


-P33 


-Q+ sin 2^12523 


-(Q2V?2 + Q|3S?2) sin 2^23 


RePai 


2Q|2COs26'i2C23 


(Q^3-2y sin 2^12^23 


RePsi 


-2 Q+ COS 2^12^23 


(Q^3-2^3)sin2^i2C23 


ReP32 


Q+2 sin 2^12 sin 2^23 


2(Q2V?2 + Q^3^?2) COS 2^23 


ImPsi 


45i2C23 


2 (523-^13) sin 26*12523 


ImPsi 


-45i2S23 


2 (^23 - ^13) sin 26*i2C23 


ImP32 





4(523C?2+5i3S?2) 






647r2 ej^ 






Pll 









P22 


(^^3 - Atsj sin 26*12 sin 26*23 






P33 


- [^^3 - Ats^ sin 26*12 sin 26*23 






RePsi 


4tec?2 + ^^3S?2)s23 






ReP3i 


4(^|3C?2+^^3^?2)C23 






ReP32 


2 ('^+3 - ^+3) sin 2^12 cos 2^23 






ImP2i 


4 &C?2 + ^23^?2) 523 






ImPsi 


4 (^r3C?2 + ^2"3S?2) C23 






ImP32 


2(S23-^r3)sin26^i2 





Table 6: Coefficients of Pfg in the RG evolution equations of the mixing 
angles 6*12, 6*^3 and 6*23, in the limit 6*13 ^ [571 fTOTj . 



B^ 
"13 



B^ 
"23 



4 {ml + mg) Jcp ± Smims{Jcp cos 20i — J'^p sin 20i) 



«^"^atm 



1 + 



4 (rrig + m|) Jcp ± 8m2ms{Jcp cos 202 — -^cp sin '■ 



(141) 
(142) 



where a = S12C12S23C23. Clearly these quantities depend on Jcp, J^p i^ 
addition to the masses and Majorana phases and hence are basis-dependent. 
In the Vj basis, all the the quantities appearing in the evolution equations 
(11391) - (I142P have finite well-defined limits for ^13 — )■ and so will be 6*13 at 
any energy scale. 

Table [7] shows the generic enhancement and suppression factors [Ff] for 
the evolution of the mixing angles, which is useful to estimate the RG evo- 
lution effects on the angles, when the active neutrinos are quasi-degenerate 
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^12 


^23 




d. 


n.h. 


i.h. 


d. 


n.h. 


i.h. 


Pll 
P22 

P33 

RePsi 

RePgi 

ReP32 
ImPai 
ImPsi 
ImP32 


Am|, 

Arrtg 
Amg 

m2 
Arrtg 

m-^ 
Arrig 

m2 
Am^ 

0(^13) 


0(^13) 
0(^13) 


0(^13) 


0(^13) 

m2 


0(^13) 
1 
1 

1 


1 
1 

1 

0(^13) 
0(^13) 


2 


■m2 




2 


AmLm 

2 


m2 


2 


Am^,„ 



Pll 

P22 

P33 

RePsi 
ReP3i 
ReP32 
ImPai 
ImPai 
ImPaa 



^13 



0{9^ 



13J 



Am 



Am 



atm 

m n 



ns 



13 



m^ a 

A^^l=^ 

m n 
A^^13 

m ri 
A^^13 

rrt /I 

A^^l=^ 

m^ /I 
A^^13 

m a 



A"^ 



13 



n.h. 



(^(^?3) 
V^^13 

a/C^13 

^13 

^13 
a/C^13 

^13 

^13 
a/C^13 



i.h. 



C^(^?3) 



'13 



'13 



(^(^?3) 
^13 



'13 



o{ei,) 



Table 7: Generic enhancement and suppression factors for the evolution of the 
angles, yielding an estimate of the size of the RG effect |S7]. The table entries 
correspond to the terms in the mixing parameter RG evolution equations 
with the coefficient given by the first column. A '1' indicates that there is no 
generic enhancement or suppression, 'd.' stands for a degenerate neutrino 
mass spectrum, i.e. Am^^^^ ^ m\ 
normally hierarchical spectrum, i.e. 



rrin 



m\ ~ 



m^. 'n.h.' denotes a 
mi <^ 1712 ^ ^^3; and 'i.h.' means an 



inverted hierarchy, i.e. m^, <^ nii < m2. 
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(Am^tjn <^ mf ~ ml ~ m| ~ m^ and this case is denoted by 'd.'), or have 
normal mass hierarchy {mi ^ m,2 ^ m,^ and denoted by 'n.h.') or inverted 
mass hierarchy (7713 ^ ?7ii < m,2 and denoted by 'i.h.'). From Table [7] we see 
that all terms in Ou are enlarged by ml^/Am!^ for quasi-degenerate masses. 
Thus, there will be large RG effects, if the different terms do not cancel each 
other. The term involving ImP32 is an exception, because its leading order 
is proportional to ^13, so that it only plays a role in special cases. Also the 
terms involving ImP2i and ImP3i will have small contributions for small 
values of (20i — 202)- In the case of a strong normal hierarchy, there is no 
enhancement. For an inverted hierarchy, where the evolution is generically 
enhanced by (~^, because the masses rrii and m2 are almost degenerate. 

Both for 6*23 and Of^,, the evolution does not depend on Pn for Oi^ = 0. 
For these two mixing angles, the enhancement and suppression factors are 
similar (for Of^ evolution, there is always an extra factor of ^13 compared 
to 6*23, as expected). The terms proportional to the other Pfg are enhanced 
by m^/Am^^j^ in the degenerate case, so that effects are expected to be 
significant, but smaller than 612 running. For both hierarchical spectra, the 
running is slow, as can be seen from Table [71 In case of diagonal P (or with 
P32 as the only non-zero off-diagonal entry) and inverted hierarchy, there will 
be no running for Of^ if ^13 = 0. However, this is no longer true if P21 or P31 
is non-zero. 

Thus in the evolution equations of the mixing angles, the generic char- 
acteristics of the terms which are proportional to the diagonal elements of 
P in the high energy theory is the same as those in the low energy effective 
theory, as already discussed in Eqs. (172l) - fl74l) in Sec I3.1.1[ 

If the diagonal elements are equal, their contributions to the RG evolution 
equations cancel exactly. This follows from the fact that the mixing angles 
do not change under RG evolution, if P is the identity matrix and thus does 
not distinguish between the flavors. As can be seen in the next few sections, 
this statement holds also for the RG evolution of the CP phases. It provides 
a consistency check for the results. Interesting new effects occur for non- 
zero off-diagonal elements in P. Some of their coefficients in the evolution 
equations do not vanish for vanishing mixings, e.g. the coefficient of P21 in 
612 in Table [HI and thus non-zero mixing angles are generated radiatively. 
This is in striking contrast to the region below the see-saw scale, as can be 
checked from Eqs. (I72|) - (l74l) . 
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647r^ Jcp/a 


647r^ J'^p/a 


Pll 








-P22 


-4a^; 


2a(^o~ - 2^-) 


-P33 


4a^7 


-2a(^o- - 2^-) 


ReP2i 


4^23^+ 


2^23(^0^ + 2^+) 


RePai 


4C236^+ 


2c23{g^ + 2g^) 


ReP32 


-2 sin 2^12 cos 2^23 Q^ 


sin2^i2Cos2^23(^(r-2^c") 


ImPsi 


2s2s{g^-2g+) 


4523^^+ 


ImPsi 


2c23{g^-2g+) 


4C236^+ 


ImP32 


sin2^i2(^o"+2^-) 


-2 sin 2^12^7 



Table 8: Coefficients of Pfg in the RG evolution equations of the Jarlskog 
invariant Jqp, the quantity J^p = Jcpcot(5, in the limit 6*13 — )• 0. The 
convention used here is a = S12C12S23C23, and Jqp = (a/2)si3C^3sin5 |1U7] . 

4.2.2 Evolution of Jcp, J'cp 

The coefficients for the RG evolution of Jcp and J'^p are presented in Table [HJ 
where the quantities g^^ ^ are given by 



Gt 
Gt 



■minis sin! 



± 



Aml^il + C) 



ml + mg 

A"^atm(l + C) ' 
)l 777,27723 sin 202 



777i7773 COS 20i 77727773 COS 2 

A"^atm(l + C) ^ 



A"^atm 



(143) 
(144) 

(145) 



Thus the quantities defined in Eqs. fll43p -( 1T^5l) are functions of masses and 
Majorana phases and hence are well-defined at all energies and at every point 
in the parameter space. Thus Table [8] shows that the running of Jcp, Jcp 
does not depend on themselves and hence independent of the Dirac CP phase 
6 upto 0{9l^). It also shows that if P is identity (or proportional to identity), 
there will be no RG evolution, as expected. 

From the generic enhancement and suppression factors for the RG evolu- 
tion of Jcp and Jcp given in Table M it can be seen that for degenerate light 
neutrino masses the coefficients are enlarged by the factor 7r7^/A7r7^^j^, for all 
Pfg except Pii. The leading contribution from Pn comes only at O^Oi^) and 
is also independent of the mass ordering of the neutrinos. For Jcp, the con- 
tributions from the other two diagonal elements P22 and P33 are suppressed 
by ^/( for normal hierarchy, while for inverted hierarchy the leading contri- 
bution is only at O^Oi^). For Jcp the evolution is suppressed by ( in both 
the cases. 
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Jcp/a 


Jcp/a 




d. 


n.h. 


i.h. 


d. 


n.h. 


i.h. 


Pll 

P22 

-P33 

RePsi 

ReP3i 

ReP32 

ImPai 
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ImP32 


0(^13) 

m2 


0(^13) 

1 
1 

c 


0(^13) 
0(^13) 

0(^13) 

0(^13) 

0(^13) 

0(^13) 

1 

1 

c 


0(613) 

2 


0(^13) 

c 
c 

1 
1 

c 


0(^13) 

c 
c 

1 
1 

c 
0(^13) 

0{9i3) 

0(613) 




rn 


m 






m 




rn 


m 






m 


Am2 
m 




Am^,^ 


Am^,^ 



Table 9: Generic enhancement and suppression factors for the evolution of 
Jcp and Jqp, yielding an estimate of the size of the RG effect. 



From Table [8] it can be seen that even if ^13 is zero to start with so that 
Jcp = Jqp = at the high scale, it can be generated radiatively. This is true 
even with a diagonal P, if P is not proportional to identity. This happens in 
the low energy effective theory also, as can be seen from Eqs. fll02|) -f lT03|l . 
However, if the diagonal elements are equal, their contributions cancel exactly 
and there will be no running at all. Table [8] also suggests that for non-zero 
off-diagonal elements of P, Jcp and J^p can be generated radiatively even 
when all the mixing angles are zero at the high scale and this is very different 
from what is expected in the effective theory below the seesaw scale. 



4.2.3 Evolution of Majorana phases 

The expressions for the running of the Majorana phases are the same in Vj 
and Vs- Table [TO] shows the running of the difference between the Majorana 
phases |0i - 02| FTIITUT]. 

As can be seen from Table [TTl the generic enhancement factors for the 
RG evolution of (cf)! — (j)2) are very similar to those for the running of 612, 
for degenerate light neutrino masses or for an inverted hierarchy. For normal 
hierarchy, there is no running if P is real, upto the zeroth order of ^13, 
which implies that each of the Majorana phases runs by equal amount. The 
running of individual Majorana phases is discussed in [57]. The running of 
the Majorana phases is also important to understand the evolution of the 
mixing angles, since all the quantities defined in Eqs. (I137p -( 1T42|) depend 
on the Majorana phases. RG evolution of the Majorana phases controls the 
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327r2(0i-02) 


Pn 


-45i2 COS 2^12 


P22 


4i5i2C23 COS 2^12 


-P33 


45i2S23 COS 26^12 


RePai 


— 81S12C23 COS 26'i2 cot 26*12 


RePsi 


85i2S23 COS 26'i2 cot 26'i2 


ReP32 


-4>Si2 COS 2^12 sin 2^23 


ImPsi 


-4Qr2C23 cot 2^12 


ImPgi 


4Qr2'S23 cot 26*12 


IniP32 






Table 10: Coefficients of Pfg in the RG evolution equations of the Majorana 
phase difference (0i — 02), in the limit 6*13 — )■ [571 HOT] . 





4>1 - 02 




d. 


n.h. 


i.h. 


Pll 

-P22 

P33 

ReP2i 
RePai 
ReP32 
ImPsi 
ImP3i 
ImP32 


Ami 

m2 
Am^ 

m 
Am'p 

m2 
Ami 

m2 
Am^ 

m 
Ami 

m2 
Ami 

m2 
Ami 


0(^13) 
0(^13) 

0(^13) 
0(^13) 
0(^13) 

1 

1 
0(^13) 


0(^13) 



Table 11: Generic enhancement and suppression factors for the evolution of 
the difference of Majorana phases (0i — ^2) [SZj- 

running of Jqp, J'cp also. 

4.2.4 Evolution of light neutrino masses 

Table [12] shows the RG evolution of the light neutrino masses. As can be 
seen, the coefficients are independent of Jcp, J'cp and hence the expressions 
remain the same in the basis P5. 

As is clear from Table [T^ the evolution of rrii is proportional to itself. This 
is a general characteristic of the running of the mass eigenvalues at all energy 
scales. As a consequence, the mass eigenvalues can never run from a finite 
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IGn'^nii/mi 


167r^m2/r«2 


IQn'^rn-i/ms 


au 


1 


1 


1 


Pll 
P22 
P33 
ReP2i 


2c?2 

0„2 2 
^■^12 '-23 
9,2 „2 
^"^12"^ 23 

-2 sin 2^i2C23 


^"^12 

0„2 2 
^'-12'-23 

2Cl2"523 

2 sin 26'i2C23 




9s2 

^"^23 

9r2 
^'-23 




ReP3i 
ReP32 
ImP2i 


2sin26'i2S23 

-2sin2^23S?2 



-2 sin 2^12^23 

-2 sin 2^23C?2 






2 sin 2^23 



ImP3i 











ImP32 












Table 12: Coefficients of Pjg in the RG evolution equations of the neutrino 
masses m, {i = 1, 2, 3}, in the limit ^13 — >^ [571 1107] ■ 



value to zero or vice versa. However, this conclusion is very specific to the 
1-loop running of the masses, and breaks down when the 2-loop contributions 
are taken into account [59]. 

As already discussed in Sec I3.1.2[ below the see-saw scales, the evolution 
of the mass eigenvalues is, to a good approximation, described by a universal 
scaling caused by the flavor-independent part of the RG evolution equations 
proportional to a,^. This flavor- independent term becomes smaller at high 
energies. Therefore, the flavor- dependent terms play a more important role 
above the see-saw scales. The importance of the flavor-dependent part in- 
creases if entries of Yx become of order one. Thus between and above the 
see-saw scales, the running may become strongly influenced by the Yukawa 
couplings of the heavy fields. 

RG evolution of Am^ and Aml^^, the quantities important for neutrino 
oscillations, can be obtained using Table [T2|, and is also discussed in |5j 



4.2.5 Contribution from Up 



(n) 



As already stated, we choose to work in the basis in which Mx is diagonal. 

(n) (n) 

Hence from the Eqs. f l20|) and fl5T]) we get that In and Ys will have non- 
zero off-diagonal components. So even if one starts with diagonal Yg [i.e. 
Yf. = Diaglne, y^jHr)) at the high scale, non-zero off-diagonal elements of Yg 

(n) (n) 



will be generated through Eqs. fllOSp and (11 101) since Y^Yx is not diagonal. 
Thus the contribution from Ue to L^pmns, as given in Eq. (167|) . will be finite 
and there will be finite contribution to the running of masses and mixing 



above and between the thresholds through F and 



ap 



Since ap is flavor 
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167r2 9^^ 


1671^91,"' 


16n^9^i 


Fn 











F22 











F33 
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-C23 


-4s23Jcp/« 





ReFai 


S23 


-4c23^Cp/« 





ReF32 
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IinF2i 





-4s23Jcp/a 





ImFsi 





-4c23Jcp/a 





ImFgs 












Table 13: Coefficients of Fjg in the RG evolution equations of all the an- 
gles {9i2, 9l^, ^23); in the limit 6'i3 — )■ 0. The convention used here is 
a = S12C12S23C23, and Jqp = (a/2) 3130^3 sin 5. We neglect ye and y^ com- 
pared to yr, and take vanishing flavor phases |107j . 



diagonal, it will contribute to the running of ye, y^ and yr, while off-diagonal 
components of F will contribute additional terms in the /3-functions of angles 
and phases. To evaluate the contributions from the off-diagonal components 
of Ye, we consider the evolution of Ue as [SZ] 



dUe 

dt 



UeX, 



(146) 



where t = ln(yu/GeV)/167r^ and X is an anti-Hermitian matrix which can be 
determined from Eqs. (166|) and (llOSp to have the form [57] 



levr^X. 



y] + yf 
yl - yf 



{UlFUe)^^ [1^3) 



(147) 



where y\ = ye and so on. The diagonal parts of X, which only influence the 
evolution of the unphysical phases, remain undetermined. Using Eqs. 
and 0146p . one can write from Eq. (I^7|) 



dU^ 



PMNS 



dt 



U, 



PMNS 



T + X^U] 



PMNS • 



(148) 



Using the expressions for T and X from Eqs. fl69|l . ( 170|l and (11471) . one gets 
the coupled equations for the angles and phases from Eq. (I148p . As suggested 
by Eq. (I148p . the contributions from the ffist term are already tabulated in 
Tables El El [lOl and [T21 and the contribution from the second term is the 
additional contribution because of the off-diagonal entries in Y^Ye generated 
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167r2 Jc^^ 


167r2 ^f-^ 


167r2 02^^ 
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F33 














ReFai 





— 523^/2 
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-C23a/2 








ReF32 














IniF2i 


-S23a/2 





C23C12/S12 


-C23 -512/012 


ImFai 


-C23a/2 





-S23C12/S12 


■S23S12/C12 


ImFaa 








-1/(C23S23) 


-1/(C23S23) 



Table 14: Coefficients of Ffg in the KG evolution equations of Jqp, J'qp and 
the Majorana phases (pi in the limit 6*13 — )■ |107j . 

in course of RG evolution. These additional terms in the /3-functions of 
angles and phases are tabulated in Tables [13] and [HI respectively. These 
contributions will just get added to the Pfg contribution for the evolution of 
the quantities given in Tables [HI M and [TDl Note that the Ffg coefficients are 
< C(l), whereas the Pfg coefficients are > O^m^/Aml^^). Since the running 
is significant only when rnf 3> Aml^^, in almost all the region of interest Pfg 
contributions dominate over the Ffg contribution^ 

In Type-II seesaw, we consider only one triplet Higgs, and so Ma is 
a number and hence from Eq. (1321) we see that Ya rnay be chosen to be 
diagonal in general. Thus if Y^, is chosen to be diagonal at high energy, it 
will remain so and the same procedure as in Sec 13.11 can be followed with the 
change that here the running of Q is to be considered instead of k. 

Note that the analytical expressions obtained in Eq. (11371) onwards, and 
those given in the tables, are valid only in the two extreme regions fi > M3 
and fi < Ml. For the intermediate energy scales, ni,^ will receive contributions 

(n) (") 

from both k and Q. In the SM these two quantities have non-identical 
evolutions, as seen from Eqs. fll35p and (11171) . and therefore the net evolution 
of Ye and laii, is rather complicated and needs numerical studies. 

Quantitative studies have been made to show that the threshold effects 
may have dramatic consequences and can make many high energy neutrino 
mass models compatible with the current oscillation data at low energy which 
would have been excluded otherwise and vice versa. To illustrate the fact, 
bimaximal mixing scenario (6'i2 = 6'23 = 7r/4, ^13 = 0) |108j is not allowed by 



^ If the running of the mixing angles 612 and 6*23 is large to make the angles close to 
zero or 7r/2 at some energy scale, the coefficients of some of the Im_Fyg may become large 
for (pi, as can be seen from Table [Hi 
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the current oscillation experiment data, as can be seen from Table [TJ But it 
is possible to make this symmetry allowed at the high energy when threshold 
effects are taken into account, in case of both Type-I |105] and Type-Ill |107j 

seesaw. 

5 Conclusions 

In the framework of the standard model (SM) of particle physics, neutrinos 
are massless at the tree level as well as at loop level. Hence one has to extend 
the SM in order to explain the tiny active neutrino masses observed exper- 
imentally. The most favored mechanisms to generate such small neutrino 
masses are the seesaw mechanisms, in which small active neutrino masses 
are generated at some high energy scale. 

All these models predict the light neutrino masses and the mixing pa- 
rameters at some high energy which corresponds, in some way, to the mass 
scale of the new fields added to the SM to generate the light neutrino masses. 
But since the experimental data are available at the laboratory energy scale, 
one needs to include the effects of renormalization group (RG) evolution. 
Unlike the quark sector where RG evolutions are quite small because of the 
hierarchical quark masses and small mixings, the effect of RG evolution on 
the neutrino masses and the mixing parameters are important. In this paper 
we reviewed the seesaw mechanisms that generate the light neutrino masses 
and the RG evolution of the neutrino masses and mixing parameters in the 
seesaw scenarios. 

The low energy effective operator to generate the light neutrino masses 
is the same for all the three types of seesaws and thus the RG evolution of 
the parameters depend only on the effective theory i.e. whether it is the 
SM or the Minimal Supersymmetric Standard Model (MSSM), which has 
been discussed in detail. The Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
parametrization of the neutrino mixing matrix is characterized by the fact 
that at ^13 = 0, the Dirac CP phase S is unphysical. This leads to the singular 
behavior of 6 at ^13 = 0. However, this singularity is unphysical, since all the 
elements of the neutrino mixing matrix t/pMNS are continuous at ^13 = 0, and 
in fact the value of 6 there should be immaterial. The singularity also creeps 
in the running of 6*13, while evolution of the other parameters is well-behaved. 
However, as discussed here, it is possible to express the RG evolution of all the 
parameters as continuous first-order differential equations if one chooses the 
basis to be Vj = {rrii, 9i2, ^23, ^13, <^i, Jcp, -^cp}' instead of the conventional 
Vs = {rrii, 9ij, (j)i, 6} basis. The evolution equations in the new Vj basis and 
their approximate integrated forms have also been discussed. 
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When the RG evolution in the high energy theory is considered, one 
needs to take the effects of the heavy fields into account carefully, since the 
evolution will depend on their interactions with the other fields. Moreover, 
the heavy fields will decouple from the theory step by step at their respective 
mass scales and start contributing through the effective operators, so the 
threshold effects are to be considered and matching conditions are to be 
imposed. We considered the evolution in three seesaw scenarios and with 
the SM as the effective theory. For energy regimes higher than the mass of 
the heaviest particle and lower than the lightest one, the evolution can be 
expressed by simple analytic formulae and qualitative understanding of the 
RG evolution is possible independent of the low energy theory considered. 
However, the final evolution of any parameter will depend on the choice of 
the high energy seesaw scenario, as well as the low energy effective theory. 

RG evolution can have dramatic effects on the masses and the mixing 
parameters in the neutrino sector, especially when threshold effects are taken 
into account. These effects can make many high energy neutrino mass models 
compatible with the current oscillation data at low energy which would have 
been excluded otherwise and vice verse. Since precision data is expected 
from the upcoming neutrino experiments, it is important to consider the RG 
evolution effects while talking about the neutrino mass models, and it then 
be possible exclude different classes of models and to gather knowledge about 
the possible high scale symmetries. 
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Appendix: Diagonalization of neutrino mass 
matrix 

To check the diagonalization procedure, let us consider the case when there 
are arbitrary 'n' number of heavy right-handed neutrinos and three active 
neutrino species. Then the neutrino mass matrix A^,^ is a n x n matrix given 
by Eq. (^ as 

M.= ( \ ^M , (A.l) 
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where mn is a 3 x n matrix, and Mn is n x n. Let us now consider the 
unitary transformation 



A B \^ f TmD\fA BV /mi 



C D J \m.l M:^ J \C D J \ ma 
tliat diagonahze the mass matrix. From Eq. flA.2p one gets 



(A.2) 



mi = AVdC* + C^jnlA* + C^M^C* , (A.3) 

m2 = B^joidD* + D^joilB* + D^M^D* . (A.4) 

The relation in Eq. flA.4p gives 

DMi2D^ = Az) + (1 - e)MN(l - e^) , (A.5) 

where A^ = D^WdII — e"^) + (1 — e)m^-B*D"^. Here we have used the fact 
that Mn ^ m^i and defined DD^^ = 1 — e, where e = (m^Mj;^^) , ra to be 
determined. The unitarity condition gives B'^B + D'^D = 1, which imphes 
B = (9(e^/^). The same holds for C. In a similar way one has from unitarity 
AA'f = 1 - e. Eq. ([Qj) also gives 

A^jcnoD* + C^m^B* + C^M^D* = , (A.6) 

which reduces to 

AVd + C^Mn = C^TnlB*{D*y^ 

= 0{e'/')inlO{e'/')D^{l + e''). (A.7) 

Thus 

Ct = 0(e) - AVdMn^ , (A.8) 

which in turn shows that e ~ 0{{j[n£)/M.-^) ). Hence B ^ C ^ (^(m/j/MN) 
and Dm2-D"^ = MN + 0(in|,/MN), which is the same as that given in Eq. f lT2|) 
in Sec[2J From Eq. (lA.Sp . keeping terms upto C(m£)/MN), one gets 

AjaiA'^ = -midM^^jo^ , (A.9) 

and this is nothing but the seesaw relation, quoted in Eq. ( ITT]) . Related 
discussions can also be found in |1U9] . 



47 



References 

[1] Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 
81, 1562 (1998). 

[2] Ahn, M. H. et al. Phys. Rev. D74, 072003 (2006). 

[3] Michael, D. G. et al. Phys. Rev. Lett. 97, 191801 (2006). 

[4] Schwetz, T., Tortola, M., and Valle, J. W. F. New J. Phys. 10, 113011 
(2008). 

[5] Apollonio, M. et al. Eur. Phys. J. C27, 331-374 (2003). 

[6] F. Dydak et al, Phys. Lett. B 134, 281 (1984). 

[7] P. Astier et al. [NOMAD Collaboration], Phys. Lett. B 570, 19 (2003) 
|a?Xw:hep-ex/0306037 | . 

[8] Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 87, 071301 
(2001) | arXiv:nucl-ex/0106015 |. 

[9] Fukuda, S. et al. Phys. Rev. Lett. 86, 5656-5660 (2001). 

[10] Bahcall, J. N., Gonzalez-Garcia, M. C, and Pena-Garay, C. JHEP 07, 
054 (2002). 

[11] T. Araki et al. [KamLAND Collaboration], Phys. Rev. Lett. 94, 081801 
(2005). 

[12] B. Pontecorvo, Sov. Phys. JETP 7, 172 (1958) [Zh. Eksp. Teor. Fiz. 
34, 247 (1957)]. 

[13] B. Pontecorvo, Sov. Phys. JETP 26, 984 (1968) [Zh. Eksp. Teor. Fiz. 
53, 1717 (1967)]. 

[14] V. N. Gribov and B. Pontecorvo, Phys. Lett. B 28, 493 (1969). 

[15] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28, 870 
(1962). 

[16] G. L. Fogli, E. Lisi, A. Marrone, A. Palazzo and A. M. Rotunno, Phys. 
Rev. Lett. 101, 141801 (2008) |arXiv: 0806. 26491 [hep-ph]]: 

[17] A. Bandyopadhyay, S. Choubey, S. Goswami, S. T. Petcov and 
D. P. Rov. larXiv:0g04r4857 [hep-ph]. 

48 



[18] M. Maltoni and T. Schwetz, larXiv:0812.3T6T] [hep-ph]. 

[19] G. L. Fogli, E. Lisi, A. Marrone, A. Palazzo and A. M. Rotunno, Phys. 
Rev. Lett. 101, 141801 (2008) [arXiv:0806.26 49 [hep-ph]]. 

[20] Kraus, C. et ah Eur. Phys. J. C40, 447-468 (2005). 

[21] R. G. H. Robertson [KATRIN Collaboration], J. Phys. Conf. Ser. 120, 
052028 (2008). 

[22] Klapdor-Kleingrothaus, H. V. et al. Eur. Phys. J. A12, 147-154 
(2001). 

[23] Komatsu, E. et al. Astrophys. J. Suppl. 180, 330-376 (2009). 

[24] S. Hannestad, Phys. Rev. Lett. 95, 221301 (2005) 
[arXiv:astro-ph/0505551J . 

[25] R. N. Mohapatra and J. W. F. Valle, Phys. Rev. D 34, 1642 (1986). 

[26] A. Zee, Phys. Lett. B 93, 389 (1980) [Erratum-ibid. B 95, 461 (1980)]; 

[27] R. N. Mohapatra and P. B. Pal, "Massive neutrinos in Physics and 
Astrophysics" (Third Edi.), World Scientific, 2004. 

[28] A. Zee, Phys. Lett. B 161, 141 (1985); A. Zee, Nucl. Phys. B 264, 99 
(1986); K. S. Babu, Phys. Lett. B 203, 132 (1988). 

[29] K. S. Babu, S. Nandi and Zurab Tavartkiladze, larXiv: 0905.27X01 [hep- 
ph]. 

[30] S. Weinberg, Phys. Rev. Lett. 43, 1566 (1979). 

[31] Minkowski, P. Phys. Lett. B67, 421 (1977). 

[32] Yanagida, T. Proceedings of the Workshop on the Unified Theory and 
the Baryon Number in the Universe. KEK, Tsukuba, Japan, (1979). 

[33] Gell-Mann, M, R. P. and Slansky, R. Complex spinors and unified 
theories. Supergravity. North Holland, Amsterdam, (1979). 

[34] Glashow, S. L. The future of elementary particle physics. Proceedings 
of the 1979 Cargese Summer Institute on Quarks and Leptons. Plenum 
Press, New York, (1980). 

[35] Mohapatra, R. N. and Senjanovic, G. Phys. Rev. Lett. 44, 912 (1980). 

49 



[36] J. Kersten, Diploma Thesis (2001). 

[37] Magg, M. and Wetterich, C. Phys. Lett. B94, 61 (1980). 

[38] Lazarides, G., Shafi, Q., and Wetterich, C. Nucl. Phys. B181, 287 
(1981). 

[39] Foot, R., Lew, H., He, X. G., and Joshi, G. C. Z. Phys. C44, 441 
(1989). 

[40] Ma, E. and Roy, D. R Nucl. Phys. B644, 290-302 (2002). 

[41] M. A. Schmidt, Phys. Rev. D 76, 073010 (2007) [ arXiv:0705.384T] [hep- 

ph]]. 

[42] M. A. Schmidt, Diploma Thesis (2004). 

[43] W. Chao and H. Zhang, Phys. Rev. D 75, 033003 (2007) 
|arXiv:hep-ph/0611 323|. 

[44] I. Gogoladze, N. Okada and Q. Shafi, Phys. Rev. D 78, 085005 (2008) 
[ arXiv:0802.3257l [hep-ph] ] . 

[45] Bajc, B. and Senjanovic, G. JHEP 08, 014 (2007). 

[46] Bajc, B., Nemevsek, M., and Senjanovic, G. Phys. Rev. D76, 055011 
(2007). 

[47] Abada, A., Biggio, C, Bonnet, F., Gavela, M. B., and Hambye, T. 
JHEP 12, 061 (2007). 

[48] Franceschini, R., Hambye, T., and Strumia, A. Phys. Rev. D78, 033002 
(2008). 

[49] del Aguila, F. and Aguilar-Saavedra, J. A. Nucl. Phys. B813, 22-90 
(2009). 

[50] Abada, A., Biggio, C., Bonnet, F., Gavela, M. B., and Hambye, T. 
Phys. Rev. D78, 033007 (2008). 

[51] Ma, E. and Suematsu, D. Mod. Phys. Lett. A24, 583-589 (2009). 

[52] Fileviez Perez, P. Phys. Rev. D76, 071701 (2007). 

[53] Dorsner, I. and Fileviez Perez, P. JHEP 06, 029 (2007). 



50 



[54] Mohapatra, R. N., Okada, N., and Yu, H.-B. Phys. Rev. D78, 075011 
(2008). 

[55] M. Malinsky, T. Ohlsson and H. Zhang, Phys. Rev. D 79, 073009 (2009) 
[ arXiv:0903.1 96T] [hep-ph] ] . 

[56] M. Mahnsky, T. Ohlsson, Z. Z. Xing and H. Zhang, Phys. Lett. B 679, 
242 (2009) [ arXiv:0905.2889 [hep-ph]]. 

[57] S. Antusch, J. Kersten, M. Lindner, M. Ratz and M. A. Schmidt, JHEP 
0503, 024 (2005) | arXiv:hep-ph/0501272j ; 

[58] S. Antusch, M. Drees, J. Kersten, M. Lindner and M. Ratz, Phys. Lett. 
B 525, 130 (2002) |arXiv:hep-ph70110366| . 

[59] S. Davidson, G. Isidori and A. Strumia, Phys. Lett. B 646, 100 (2007) 
|arXiv:hep-ph/0611 389] . 

[60] M. E. Peskin and D. V. Schroeder, "An Introduction to Quantum Field 
Theory," (Addison- Wesley, Reading, Massachusetts, 1997). 

[61] S. Antusch, J. Kersten, M. Lindner and M. Ratz, Nucl. Phys. B 674, 
401 (2003) | arXiv:hep-ph/0305273| . 

[62] A. Dighe, S. Goswami and S. Rav. larXiv:0810.5680l [hep-ph]. 

[63] A. Dighe, S. Goswami and P. Roy, Phys. Rev. D 73, 071301 (2006) 
tarXiv:hep-ph/0602 062| . 

[64] A. Digh e, S. Goswami an d W. Rodejohann, Phys. Rev. D 75, 073023 
(2007) |arX iv:hep-ph/0612328l . 

[65] A. Dighe, S. Goswami and P. Roy, Phys. Rev. D 76, 096005 (2007) 
jarXiv:0704.3735l [hep-ph] ] . 

[66] J . R. Ellis and S . Lola, Phys. Lett. B 458, 310 (1999) 
I arXiv:hep-ph/9904'279 | . 

[67] Chankowski, P. H., Krolikowski, W., and Pokorski, S. Phys. Lett. 
B473, 109-117 (2000). 



[68] P. H. Chankowski and Z. Pluciennik, Phys. Lett. B 316, 312 (1993) 



^arXiv:hep-ph/9306333 1 ; P. H. Chankowski and S. Pokorski, Int. J. 



Mod. Phys. A 17, 575 (2002) [arXiv:hep-ph/01 10249]. 



51 



[69] K. S. Babu, C. N. Leung and J. T. Pantaleone, Phys. Lett. B 319, 191 
(1993) |arXiv:hep-ph79 309223|. 

[70] S. Antusch, M. Drees, J. Kersten, M. Lindner and M. Ratz, Phys. Lett. 
B 519, 238 (2001) |arXiv:hep-ph7010800 5|. 

[71] T. Fukuyama and N. Okada, JHEP 0211, Oil (2002) 
|arXiv:hep-ph/0205066] . 



[72] M. Tanimoto, Phys. Lett. B 360, 41 (1995) | arXiv:hep-ph/9508247| . 

[73] N. Haba, N. Okamura and M. Sugiura, Prog. Theor. Phys. 103, 367 
(2000) |arXiv:hep-ph/9810471| . 

[74] K. R. S. Balaji, A. S. Dighe, R. N. Mohapatra and M. K. Parida, Phys. 
Rev. Lett. 84, 5034 (2000) |arXiv:hep-ph/0001310| ; K. R. S. Balaji, 
A. S. Dighe , R. N. Mohapatra and M. K. Parida, Phys. Lett. B 481, 
33 (2000) |a rXiv:hep-ph/0002177j . 

[75] K. R. S. Balaji, R. N. Mohapatra, M. K. Parida and E. A. Paschos, 
Phys. Rev. D 63, 113002 (2001) |arXiv:hep-ph/0011263i . 

[76] R. N. Mohapatra, M. K. Parida and G. Rajasekaran, Phys. Rev. D 69, 
053007 (2004) |arXiv:hep-p h/0301234]. 

[77] S. K. Agarwalla, M. K. Parida, R. N. Mohapatra and G. Rajasekaran, 
Phys. Rev. D 75, 033007 (2007) |arXiv:hep-ph/0611225l . 



78] F. Vissani, |arXiv:hep-ph/9708483[ 



[79] G. C Branco, M. N. Rebelo and J. L Silva-Marcos, Phys. Rev. Lett. 



82, 683 (1999) HarXiv:h ep-ph/ 9810328j . 



[80] J. A. Casas, J. R. Espinosa, A. Ibarra and L Navarro, Nucl. Phys. B 
556, 3 (1999) |arXiv:hep-ph/9904395l . 

J. A. Casas, J. R. Espinosa, A. Ibarra and I. Navarro, Nucl. Phys. B 
569, 82 (2000) |arXiv:hep-ph/9905381|. 



[82] N. Haba, Y. Matsui, N. Okamura and M. Sugiura, Prog. Theor. Phys. 
103, 145 (2000) |arXiv:hep-ph/9908429| . 

[83] R. Adhikari , E. Ma and G. Rajasekaran, Phys. Lett. B 486, 134 (2000) 
[arXiv:hep-ph/0004197 |. 



52 



[84] A. S. Joshipura, S. D. Rindani and N. N. Singh, Nucl. Phys. B 660, 
362 (2003) | arXiv:hep-ph/0211378| . 

[85] A. S. Joshi pura and S. Mohanty, Phys. Rev. D 67, 091302 (2003) 
|a rXiv:hep-ph/0302181 ] . 

[86] Z. Z. Xing and H. Zhang, Comniun. Theor. Phys. 48, 525 (2007) 
| arXiv:hep-ph/0601106] . 

S. T. Petcov, T. Shindou and Y. Takanishi, Nucl. Phys. B 738, 219 
(2006) [ arXiv:hep-ph/0508243j . 



A. S. Joshipura, Phys. Lett. B 543, 276 (2002) [arXiv:hep- ph/0205038]. 
[89] A. S. Joshipura and S. D. Rindani, Phys. Rev. D 67, 073009 (2003) 



|arXiv:hep-ph/02 114041 . 

[90] J. W. Mei an d Z. Z. Xing, Phys. Rev. D 70, 053002 (2004) 
I arXiv:hep-ph7040408 1] . 



[91] S. Lola, Acta Phys. Polon. B 31, 1253 (2000) [ar Xiv:hep-ph/0005093| . 

[92] N. Haba and N. Okamura, Eur. Phys. J. C 14, 347 (2000) 
|arXiv:hep-ph/9906481j . 

[93] E. Ma, J. Phys. G 25, L97 (1999) |arXiv:hep-ph/9907400| . 

[94] N. Haba, Y . Matsui, N . Okamura and T. Suzuki, Phys. Lett. B 489, 
184 (2000) |arXiv:hep-ph/0 005064 |. 

[95] P. H. Chankowski, W. Krolikowski and S. Pokorski, Phys. Lett. B 473, 
109 (2000) [ arXiv:hep-ph/9910231 |. 

[96] J. T. Pantaleone, T. K. Kuo and G. H. Wu, Phys. Lett. B 520, 279 
(2001) [ arXiv:hep-ph/0108137j . 

[97] S. Lu o and Z. Z. Xing, Phys. Lett. B 637, 279 (2006) 
tarXiv:hep-ph/ 0603091) . 

[98] F. Plentin ger and W. Rodejohann, Phys. Lett. B 625, 264 (2005) 
|arXiv:hep-ph/0507 143]. 

[99] A. Dighe, S. Goswami and P. Roy, Phys. Rev. D 76, 096005 (2007) 
iarXiv:0704.3735l [hep-ph]] . 



53 



[100 
[101 
[102 
[103 
[104 
[105 
[106 

[lor; 

[108 



M. Hirsch, E. Ma, J. C. Romao, J. W. F. Valle and A. Villanova del 



Moral, Phys. Rev. D 75, 053006 (2007) |arXiv:hep-ph/0606082] . 

M. A. Schmidt and A. Y. Smirnov, Phys. Rev. D 74, 113003 (2006) 
|arXi v:hep-ph7060 7232| . 

S. Goswami, S. T. Petcov, S. Ray and W. Rodejohann, Phys. Rev. D 
80, 053013 (2009) [ arXiv:0907.286 9 [hep-ph]]. 

S. F. King and N. N. Singh, Nucl. Phys. B 591, 3 (2000) 
[arXiv:hep-ph/0006229] . 

S. Antusch, J. Kersten, M. Lindner and M. Ratz, Phys. Lett. B 538, 
87 (2002) |a rXiv:hep-ph/0203233] . 

S. Antusch, J. Kersten, M. Lindner and M. Ratz, Phys. Lett. B 544, 
1 (2002) |arXiv:hep-ph/0206078l . 

T. M iura, T. Sh indou and E. Takasugi, Phys. Rev. D 68, 093009 (2003) 
|arXi v:hep-ph/0 308109|; T. Shindou and E. Takasugi, Phys. Rev. D 70, 
013005 (2004) |arXiv:hep-ph/0402106i . 

J. Chakrabortty, A. Dighe, S. Goswami and S. Ray, Nucl. Phys. B 820, 
116 r2009l 1^X17:0812.27761 [hep-ph]. 



F. Vissani, arXiv:hep-ph/9708483j V. D. Barger, S. Pakvasa, 



T. J. Weiler and K. Whisnant, Phys. Lett. B 437, 107 (1998) 
|arXiv:h ep-ph/9806387|; A. J. Baltz, A. S. Goldhab er and M. Gold- 
haber, Phys. Rev. Lett. 81, 5730 (1998) |arXiv:hep-ph/9806540]; 
M. Jezabek and Y. Sumino, Phys. Lett. B 440, 327 (1998) 
|arXiv:hep-ph/9807310] ; G. Altarelh and F. Feruglio, Phys. Lett. B 



439, 112 (1998) |ar Xiv:hep-ph/9807353 j; 
[109] Z. Z. Xing, Int. J. Mod. Phys. A 19, 1 (2004) |arXiv:hep-ph/0307359| . 



54 



